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ELASTICA AS A DYNAMICAL SYSTEM 


LARRY BATES, ROBIN CHHABRA AND J^DRZEJ SNIATYCKI 


Abstract. The elastica is a curve in R 3 that is stationary under variations of the 
integral of the square of the curvature. Elastica are viewed as a dynamical system 
that arises from the second order calculus of variations, and its quantization is 
discussed. 


1. Introduction 

Ever since the beginning of the calculus of variations, second order problems such 
as the classical problem of the elastica have been considered. The peculiar situation 
that distinguishes most of the interesting examples in second order problems from 
the more familiar first order theory is that they are parameter independent, and 
so the theory of such problems has a somewhat distinctive tone from that of the 
more familiar first order theory. A comprehensive review of this theory, as it was 
understood up until the 1960s, may be found in the monograph by Grasser 171 W 
By way of contrast, this paper emphasizes the relation of the variational problem 
to the geometry of the corresponding Cartan form and its interplay with symmetry, 
conservation laws, the Noether theorems and the associated first order canonical 
formalism. The problem of elastica is then re-examined in light of this discussion. 

The rationale for this paper is to give a biased view of that portion of the the¬ 
ory of second order variational problems that could reasonably be expected to be 
useful for understanding the common behaviour of several geometric functionals 
on curves. Good examples to keep in mind while reading this paper (these are the 
three main examples that motivated our study) are the elastica, the shape of a real 
Mobius band in terms of the geometry of the central geodesic If20l . and the curve 
of least friction. 

For reasons not entirely clear to us, the geometric theory of higher order vari¬ 
ational problems seems to have developed in a manner largely detached from the 
needs and concerns of concrete problems. This is a startling contrast to recent 
developments in geometric mechanics and their understanding of stability, bifurca¬ 
tion, numerical schemes, the incorporation of nonholonomic constraints, etc. The 
consequences of this are at least two-fold: first, it leads to a palpable sense of 
dreaqj when faced with trying to look up a formulation of some part of the theory 
that will cleanly explain how to compute something obvious, and second, a real 
disconnect between the theoretical insights and the actual computational methods. 
This disconnect is vividly illustrated in the problem of elastica. 

'This monograph is especially noteworthy for its comprehensive bibliography. 

2 This may be reduced to mere frustration by those less ignorant than the authors. 
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Planar elastica (the equilibrium shape of a linearly elastic thin wire) were con¬ 
sidered at least as early as 1694 by James Bernoulli!] However, it was not until 
about 1742 that Daniel Bernoulli convinced Euler to solve the problem by using 
the isoperimetric method (the old name for the calculus of variations before La¬ 
grange.) From a variational point of view, the elastica is idealized as a curve that 
minimizes the integral over its length of the square of the curvature (that is, mini¬ 
mize f k 2 ds), and is thus naturally treated as a second order problem in the calculus 
of variations. Exhaustive results were then published by Euler in 1744 0 - Since 
Euler’s results were so comprehensive, it is not surprising that the study of elas¬ 
tica remained somewhat dormant until taken up again by Max Bom in his thesis 
lf2l . More recently a striking result was obtained in 1984 by Langer and Singer 
fl2l when they demonstrated the existence of closed elastica that were torus knots. 
Their proof was noteworthy because they eschewed the usual variational machinery 
and employed clever ad hoc geometric arguments such as an adapted cylindrical 
coordinate system to aid their integration. In fact, a significant motivation for this 
paper was to see to what extent their results could be understood by a more pedes¬ 
trian use of the second order calculus of variations that looked more like just ‘turn¬ 
ing the crank’ on the variational machine, and thus had the comfort of familiarity 
of technique. 

Some features of the elastica problem instantly spring to mind in the modem 
geometrically oriented reader. The first is that the problem is manifestly invariant 
by the action of the Euclidean group. The second is that it would be very nice to 
have a theory that explained how to reduce the symmetry using the concomitant 
conservation laws that Emmy Noether taught us are in the problem, and then wind 
up with some form of reduced Euler-Lagrange equations. Assuming we can solve 
these reduced equations, and hence know the curvature and torsion of the elastic 
curve, we would expect a good theory to show us methods to determine the shape 
of our curve that go beyond a trite referral to the fundamental theorem of curves 
stating that the curvature and torsion of the curve determine it up to a Euclidean 
motion. Given all this, what we actually find when we look at the published work 
on elastica (such as lfl2l or 0) is that it proceeds somewhat differently. In par¬ 
ticular, almost none of the actual computations seem to follow any method that 
resembled the current theory. There are good reasons for this, and it is not due 
to ignorance of those geometers but a reflection that the theory at that time was 
presented in such a way as to simply be unhelpful, and unable to easily identify 
the geometric meaning of some of their calculations. This is the best explanation 
we have of the situation at the time and why it was still necessary a decade after 
lIlZl appeared for Foltinek (see 0) to write a paper demonstrating the integration 
constants in elastica in terms of the conserved Noetherian momenta. 

The plan of this paper is to first discuss the Euler-Lagrange equations. This 
section will serve to fix notation and some basic notions. Next, we discuss the 
Cartan form as the geometrization of the variational problem using the structure 
of the jet bundle. This is followed by discussion of parametrization invariance and 


3 See the delightful discussion by Levien in ED or ED- 
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the Hamiltonian formalism. Elastica are then studied from this point of view. The 
problem is then recast as a constrained Hamiltonian system and the reparametriza- 
tion group and its reduction are examined. The paper concludes with the geometric 
quantization, and discusses the quantum representations of the groups SE(3) and 
Diff + R as well as the quantum implementation of constraints. 


2. Second order variational problems 


2.1. Variation of the action integral. A curve [foHi] —> R" : t i-» x(t) can be 
uniquely described by the corresponding section 

cr : [f 0 ,fi] -> [?oHi] xR":(h (1) 


We consider Q = R x R" as a bundle over R with typical fibre R' . For each integer 
k > 0, we denote by J k the A-th jet of sections of Q. Furthermore, interpret the 
section cr given in dTJ as a local section of Q and denote by j k cr : [f 0 , f| ] —> J k the 
k- jet extension of cr. This paper considers variational problems defined by second 
order Lagrangians. For a Fagrangian L : J 2 —> R : (f, x, x, x) i-» L(t, x, x, x), the 
corresponding action integral is 

rt i rh 

A(cr) = I (L o j 2 cr) dt = I L(t, x(t), x(t), x{t)) dt. 

JtQ Jt 0 

A variation of a section (without variation of time) cr cr + Scr : t x(t) + Sx(t) 
extends to the second jets as 

j 2 cr i-> j 2 cr + 6j 2 cr : t (x{t) + 6x(t), x{t) + 6x{t), x{t) + 6x{t)), 


where 

d d d 2 

6x{t) = — 5x(t) and 5x(t) = — dx(t) = —- 5x{t). 
dt dt dt 1 

Then, integrating by parts twice, it follows that the action varies as 


5A(cr) = 


~ f 

Jt 0 

- r 

'-'to 


6L(t, x(t), x(t), x{t)) dt 


{dL ddL d 2 dL 


\ dx 


1 dL d dL 
+ { M (l) -d,te WSx 


dL 

+ —(f) 5x 
dx 


to 


It follows from the fundamental lemma of the calculus of variations that 


Conclusion 2.1. The action integral 

nh nh 

A(cr) = (L o j 2 cr) dt = I L(t, x{t), x(t), x{t)) dt 

Jto Jto 

is stationary with respect to all variations cr cr + 6cr : t x(t) + dx(t), such 
that dx and 6x vanish on the boundary, if and only if the section cr satisfies the 
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Euler-Lagrange equations 

dL d dL d 2 dL 
dx dt dx dt 2 dx 


(2) 


Consider now the boundary terms in the variation. The partial derivative || is a 
map from J 2 to R", and 

dL I dL 

—(t) 5x = ( —(f, x(t), x(t), x(t)), 6x(t) 
dx \ dx 

where the angle bracket denotes the Euclidean scalar product in R' ! . However, 


d dL \ d dL d dL d dL ... d dL\ 

dt \ dx ) \ dt dx dx dx dx dx dx dx j 


depends on the third derivative x of the section cr, and hence, can be inter¬ 
preted as a map from J 3 to R”. Using the projection map 

Kyi : J —> J : (t, x, x, x, x) i-» (t, x, x, x), 

define Ostrogradski’s momenta by 



* dL 


Px 

n32 dx 



* dL 

d dL 

Px 

7132 dx 

dt dx' 


and interpret them as maps from J 3 to R". In the following, in order to simplify 
the notation, the pull-back sign is omitted and an overdot is used to denote the 
derivative with respect to t. This leads to the usual expressions 


dL 

Px = ~X7. j 

dx 

dL d_dL _ dL 
dJc dt dx dx 


With this notation, the variation equation is 


(3) 


r n (dL d dL d 2 dL } „ „ 

M(cr) = X W' ) -5?» < ' ) + ^aj <,) r ; ■<=*+p,«»+ ?»«!,,. (4) 

With an eye towards towards the Caiian form, it is convenient to reinterpret a 
variation as the Lie derivative with respect to a vector field. Let a variation cr 
i—> cr + 8cr : 1 1 —> x(t) + 8x(t) of cr be given by a vector field X on / 2 that is tangent 
to the fibres of the source map J 2 —> [to,t\\ : (t,x,x,x) i-» t. In other words, if 
X = X x 4^, then 

5x(t) = X x (cr(t)). 

Then the variation j 2 cr i-> j 2 cr + 8j 2 cr is given by the prolongation 

, d d d 
x 2 =X x —+X x —+X x — 
dx dx dx 
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of X to J 2 , where 

X t = —X x and X* = —X x = ^—X x . 
dt dt dt 2 

In other words, 

6x(t) = X x (j 2 cr(t)) and 6x(t) = X x (j 2 {crj). 


With this identification. 


5A{cr) 


r h (dL dL dL 

\—(t)6x+—(t)6x+—(t)6x}dt 
*-'tQ V ox OX ox, 

rh rh 

£ x i(Ldt) = X 2 _l d(Ldt) 

Jtn Jtn 


(5) 


because 


£ X 2(L dt) = X 2 _l d(L dt) + d(X 2 _| L dt), 

and the the assumption that X is tangent that to the fibres of the source map implies 
that X 2 _l Ldt = 0. Comparing equations © and © yields 



X 2 _l d(Ldt) 


f 

Jto 

f 

O to 


f dL d dL d 2 dL) t , 

(dL d dL d 2 dL\ / 

(e;-Jtdi + dfiTxl x ’‘ 1,+ ( p ' dx + Ptdx ' x i 


where (p x dx + p x dx,X l ) is the evaluation of a 1-form p x dx + p x dx on the first jet 
bundle on the first jet prolongation X 1 of X (see proposition (12.6b ). In equation © 
p x dx and p x dx arc interpreted as one-forms on J ] . 


2.2. The Cartan form. The contact forms of the second jet bundle J 2 are 
f)\ = dx - x dt and (h = dx - xdt. 

Their importance stems from the following 

Proposition 2.2. A section cr : [Iq, t\ \ —» J 2 : t h (t,x(t),x(t),x(t)) is the jet 
extension of the section of its projection po, l\ \ —> R" : 1 ^ (t, x(t)) by the source 
map J 2 —> [to, t\] : (f, x, x, x) (t, x) if and only ifcr*6\ = 0 and cr*0 2 = 0. 

Definition 2.3. The Cartan form corresponding to a Lagrangian L is the one-form 
0 on J 3 given by 

0 = Ldt + p x {dx - xdt) + p x (dx - xdt), (7) 

where p x and p x are the Ostrogradski momenta ©. 

Observe that 0 may by written in the form 

0 = p x dx + p x dx - H dt, (8) 


where 


H = p x x + p x x - L 


( 9 ) 
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is the Hamiltonian of the theory. Since 0 differs from the Lagrange form Ldt by 
terms that arc proportional to the contact forms, it follows that for any section cr 
the action A(cr) can be expressed as the integral of 0 over fcr. In other words, 


A(cr) = f ( Loj 2 cr)dt= f (j 2 erf Ldt = f (jV)*0. 

to 'Jt(\ d to 


( 10 ) 


Therefore, the Cartan form 0 may be used instead of the Lagrange form Ldt to 
describe the variational problem under consideration. Other aspects of the Cartan 
form are discussed in ifTOl . 

Definition 2.4. A Langrangian L is regular if the matrix 

d 2 L 


dxjdxi 


is non-singular. 


Theorem 2.5. Let y be a section of the source map J 3 —> [to, t\ \ projecting to a 
section cr of [to, 0 ] X R" —> [to, ?i] and let jfer be the the third jet extension of a. 


(1) Ify = ] 2 <t, then cr satisfies the Euler-Lagrange equations if and only if the 
tangent bundle of the range ofy is contained in the kernel of d®. 

(2) If the Lagrangian L is regular and the tangent bundle of the range of y is 
contained in the kernel of d 0, then y = y 3 cr and cr satisfies the the Euler- 
Lagrange equations. 

Proof. The exterior differential of 0 can be written as 

dL 

d® = p x dx A dt + — dx A dt + dp x A (dx - xdt) + dp x A (dx - xdt ), 
ox 

because 

d® = dL A dt - p x dx A dt - p x dx A dt + dp x A (dx - x dt) + dp x A (dx - x dt) 

= (L x - p x ) dx A dt + (L x - p x ) dx A dt + L x dx A dt + 

+dp x A (dx - xdt) + dp x A (dx - xdt) 

= p x dx A dt + L x dx A dt + dp x A (dx - x dt) + dp x A (dx - x dt). 

A section y : t i-» y(t) = (t, x(t), x(t), x(t)) of the source map projects to a section 
cr : t cr(t) = (t, x(tj). Moreover, 

. d dx(t) d dx(t) d dx(t) d 
Ty(d,) = — + + 


dt dt dx dt dx dt dx, 
spans the tangent bundle of the range of y, and 

Ty(d t )(p x ) - p x and Ty(d t )(p x ) - p x . 
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Hence, 


Ty(d,)-\d® = 

. dtft) . .. dLdx(t) dL 

= px —-— at - p x dx + —-- —dt - — dx + p x (ax - x at) + p x (dx - xdt) + 

dt ox dt ox 

- (dx - xdt, Ty(d t )) dp x - (dx - x dt, Ty(d t )) dp x 


I dx(t) 
P *\-* 


. , dL (dx(t) 

* U,+ -A— 


dL 

x ) dt + ( p x - — | (dx - xdt) + 

OX ) 


- (dx - xdt,Ty(d t )) dp x - (dx - xdt, Ty(d t )) dp x . 

To prove the first statement observe that if y = fcr, then 

dx(t) 


dt 

dx(t) 

dt 


- x = 0, 


-xdt = 0, 


(dx - xdt, Ty(d t )) = 0, 

(dx - xdt, Ty(d t )) = 0, 


and 


C\T \ 

Ty(d t ) _l d@ = \ p x - — j 


(dx - xdt). 


Hence, Ty(d t ) is in the kernel of dQ if and only if p x ~ 'jp. = 0. By definition of the 
Ostrogradski momenta. 


dL 


(dL 

d dL\ 

dL 

dx 

dt 


dt dx ) 

dx ’ 


Px 

and so Ty(d,) is in the kernel of d@ if and only if the section cr satisfies the Euler- 
Lagrange equations. 

To prove the second paid, suppose that Ty(d t ) is in the kernel of dQ. Then 

dL) 


( dx(t) 

Pi \dT 


. , dL (dx(t) 

hr 


x]dt + \p x - — \(dx - xdt) + (11) 

dx i 


- (dx - xdt, Ty(d t )) dp x - (dx - xdt, Ty(d t )) dp x = 0. 

Evaluating the left hand side on a vector v4», where v is an arbitrary vector in R", 
yields 

(dx - xdt, Ty(d t )) v-^- = 0, 
dx 

because only p x depends on x. Since 


dL d dL dL d 2 L . d 2 L 

dx dt dx dx dtdx 


d 2 L 


d 2 L 


- x 


- x- 


dxd'x dxd'x d'xd'x ’ 


it follows that 


dp x _ d 2 L 
dx dxd'x 


n d 2 L 

(dx - xdt, Tcr(d,)) v —— = 0 
d'xd'x 


Therefore, 
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for an arbitrary vector v. By assumption of regularity of the Lagrangian, the matrix 
(Jm?) is non-singular. Hence, 

(dx - xdt, Tcr(d t )) = 0, 


which implies that - x(t) = 0. Substituting these results into equation (fTTb 
yields 

- x\dt + (p - — - xdt) - (dx - xdt,Ty(d t )) dp* = 0. (12) 


Px 


I dx(t) 
\ dt 


Evaluating the left hand side of this equation on a vector vj|, where v is an arbitrary 
vector in R", yields 


(dx - xdt, Ty(d t )) = 0. 

ox 


Since 


dpx 

dx 


d 2 L 

d'xd'x 


d 2 L 


(dx - xdt,Ty(d t )> v—— : = 0 
oxox 

for every vector v in R'. The assumed regularity of L implies that 

(dx - xdt , Ty(d t )) = 0, 

so that - x = 0. Hence, y - fcr and equation (fl2T) reads 

8L\ 

px - — I (dx - xdt) = 0, 
which implies that cr satisfies the Euler-Lagrange equations. 

2.3. Symmetries and conservation laws. 


q.e.d. 


2.3.1. Symmetries of the Lagrange form. Consider an infinitesimal transformation 
in (to, t\) x R'' given by 

t = t + er(t, x), x l = x 1 + ef‘(t, x). 

It corresponds to a local one-parameter group of local diffeomorphisms generated 
by the vector field 

X = T—+f‘—. (13) 

dt % dx 1 

Proposition 2.6. The prolongations X 1 , X 2 and X 2 of the vector field X in equation 
(El to the jet bundles J ] , J 2 and J 3 , respectively, are 

• d : d d 

x l = T - \-L 7 + (f - XT ) , 

dt dx 1 dx 


d 


d 


d 


d 


X = t— + p—+ (f - xf)—+ (f - 2xt - xf)—, 


dt 

d 


dx 

d 


dx 

d 


dx 

d 


X = T—+f— + (f-xt)— + (f-2xT-xf)— + (p-3xt-3xf-xf)—. 

dt dx dx dx dx 
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It remains to relate the prolongations of X to the contact forms 0\ - dx - xdt 
and 82 = dx - x dt. 


Proposition 2.7. Let I = [to, ti]. For a section cr ofl x R" —> 1, 
j l cr*£ x i9i = 0 and j 2 cr *£ X 202 = 0. 


Proof. Since 


then 


i d d 

X 1 =T-+f- + (f-xr)-z7 


dt dx 


dx 


£ x \ 0 \ = X 1 _I dG\ + d{ 6 [ ,X 1 > 

= \-r--x-r\ idx-xdt). 
\ dx dx 


Therefore 


j l cr*£ x i 6 \ = j l cr*(-f dt + xf dt + df - xdf), 
= -jj dt + xf dt + f dt - xf dt, 

= 0. 

On the other hand, since 62 = dx - x dt, and 


d . d 


d 


X- = T—+ f— + (£- XT)— +Q- 2 XT - xf)— 


dt “ dx 


dx 


dx 


it follows that 


£ x 2 0 2 = X 2 _| d0 2 + d(0 2 ,X 2 ) 


(df - fdij - T{dx - xdt) - x^—, (dx - xdt) 

I d (dx dx 

+x ( — — + — x 


dx 


dx \ dt dx 


dx 

x 1, (dx - xdt)) + x ( —, (dx - xdt) 
dx 


Therefore, 


j 2 (T*£ x i02 = j 2 cr* i^df - f dtj - t(dx - xdt) - xl^—,(dx - xdt)j^ + 


+j 2 cr* i x 


d Idx dx 

— — + —x 


dx \ dt dx 


x\, (dx - xdt)) + i ( —, (dx - xdt) 


dx 


dx 


= 0. 


q.e.d. 


Let cr be section of / x R" —> /, where / = [to, ti] and 


A (a) = J"(L° fcr)dt = J' j 2 cr*(Ldt) - J" 


Ldt. 


(14) 


j-cr(I) 
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Proposition 2.8. The action integral ( 1/-/D is invariant under the one-parameter 
local group exp tX 2 of local diffeomorphisms of J 2 generated by X 2 if 

j 2 cr*(£ X 2 (Ldt)) = 0. 

Moreover, ^Afexp tX(cr))\ t= o for every section cr, if and only if £^ (Ldt) = 0. 


Proof Consider a vector field X on / x R' 1 . Denote by exp tX the local one- 
parameter group of local diffeomorphisms of 7xR" generated by X, and by exp tX 2 
the local one-parameter group of local diffeomorphisms of J 2 generated by the pro¬ 
longation X 2 of X to J 2 . Then, 

A(exp tX(cr)) - I Ldt. 

Jexp tX 2 U 2 o-(P) 

The change of variables theorem asserts that ( ( co - £ dfoo for any form to, chain 
c and a one-parameter group tp, diffeomorphisms. Hence, 


Therefore, 


f Ldt = f 

Jexp tX 2 {j 2 cr(I)) Jj-ail ) 


Ldt = I (exp tX 2 )*Ldt 


Afexp tX(cr)) = f (exp tX 2 )*Ldt. 
df-o-(T) 


Now, differentiating under the integral sign with respect to t, 


7 i /-* 

—Afexp tX(cr)) = I -^-(exp tX 2 )*Ldt = I (exp tX 2 )*£ x iLdt, 

ut jp-cr(J) dt JJ2 q-(D 


and setting t = 0, 


^-A(exp tX(a))\ t= Q = f £ x iLdt = f j 2 cr*(£ x 2 (Ldt)), 
dt J j 2 (r(I) JI 

it follows that if j 2 cr*£ x iLdt = 0, then ^Afexp ?X(cr))| f= o = 0. Moreover, 
^Afexp tX(cr ))| f=0 = 0 for every section cr, if and only if £ x iLdt = 0. q.e.d. 


Definition 2.9. A vector field X on I x R" is an infinitesimal symmetry of the 
Lagrange form Ldt if £ x i(Ldt) = 0. 

Lemma 2.10. The Lie derivative of the Lagrange form Ldt with respect to X 2 is 


/ UL-/ ■ UJ-o U1 -j UL-j \ 

£x 2 (Ldt) = It— +£ — + (f-xf)— + (f - 2xt - xf) — \dt + L dr. 

Hence, for every section cr of I X R" —> I, 

, * I dL ML dL .. dL \ 

j-cr (£ X 2 (Ldt)) = It— + f — + (f - xt)— + (f - 2 xt - xt)— + Lt\ dt. 
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Proof. The Lie derivative of the Lagrangian form L dt with respect to X 2 is 
£ X 2 (L dt) = X 2 _l dL A dt + d(X 2 _| L dt) 

, ((dL dL dL dL \ \ , 

= X _l — dt +—rdx + —dx. + —d.r A dn + t/(X“ _l Ldt) 

\\ dt dx 1 dx dx ) ) 

t ((dL dL dL \ \ 0 

= X _l —-dx H-dx H- dx A dn + d(X“ _l L dt) 

\\dx’ dx dx / 


because dt A dt = 0. Hence, 
£ X 2 (Ldt) - 


:dL . dL .. dL\ 

f ~ xt)— + (f- 2 xt - xt)— dt + 

dx' dx dx 


(dL dL dL \ 

-t —: dx + —dx + — dx + L dx + t dL 
\dx' dx dx ) 

dL :dL . dL .. dL\ 

r— + £ + ^~ XT ^ + (£ ~ 2xt ~ XT )-^ dt + LdT - 

dt dx 1 dx dx 


q.e.d. 


Lemma 2.11. (Noether identities.) The equation j 1 (r*(£ x i(Ldt)) = 0 is equivalent 
to 


dt 


dL d (dL 


dL d 


-r \Lt + hr - T - (£~tx) + — — (f-Tx) \ = 


dx dt \ dix 


dL d (dL\ df (dL 


dx dt 


(15) 


Proof. The details of this routine but lengthy calculation may be found in fl6l . 

q.e.d. 

Remark 2.12. The Noether identity is essentially the extension of the equation 
j 2 cr*(£ X 2 (Ldt)) = 0 to the fourth jet bundle. More precisely, if 714,2 : d 4 —» J 2 is 
the natural projection and X 4 is the prolongation of X to d 4 , then 

A*« 2 (% i( i dim = j v ((|-|(|) + ^(|))«-”)) + <i« 


+f(T* I \LT + I ^7 - 2- I — I I (f ~ TX) + — ~ TX) 


dt 


dL d l dL 


dx dt \ dx 


dL d 


dx dt 


An immediate corollary of the Noether identity is the following conservation 
law. 

Theorem 2.13. (First Noether theorem.) To every infinitesimal symmetry X = tJj + 
of the Lagrange form Ldt, there corresponds a conserved quantity 


dL d (dL\ 


dx dt \dx ) 


dL d 


X = Lt + — - — — (£ - TX) + — — (f - TX). 


dx dt 


( 17 ) 


That is, J? x is constant along solutions of the Euler-Lagrange equations 

dL _ d_ ldL\ £_ ldL\ _ 


dx dt\dx dt 2 \ dix 
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In other words, if cr satisfies the Euler-Lagrange equations, then fa* (X' _l 0) is 
constant. 

Example 2.14. (Conservation of linear momentum.) If the Lagrangian L does not 
depend on the coordinate x\ then X = is an infinitesimal symmetry and the 
momentum 

' a? dx' dt \ dx') 

is conserved. 


Example 2.15. (Conservation of energy) If the Lagrangian L does not depend on 
the parameter t, then X = is an infinitesimal symmetry and the energy 

H - p x x + pfx - L = -'/a 


is conserved. 


Proof 
£a_ (L dt) 


d d 

— _l d(Ldt) + d(— _l Ldt) 

dt dt 

d 

— _l dL A dt + dL 
dt 


dL dL dL 
—dx + —dx + —dx | + 
dx dx dx 


(dL , dL , dL , 
—dt + —dx + — dx + 
\ dt dx dx 




Hence, ^ = 0 implies that £ is an infinitesimal symmetry of the Lagrange form 
Ldt. q.e.d. 


Other conserved quantities will be discussed later. 

Remark 2.16. There is a vast amount of work on symmetry principles and conser¬ 
vation laws following Noether’s fundamental paper liTTI . Three works in particular 
arc noteworthy: the monographs by Logan |[T6l and Kosmann-Schwarzbach |j9l, 
and a review by Krupkova liTTI . 

2.3.2. The Cartan form approach. Recall that the Cartan form 0 is 

0 = Ldt + p x 6\ + p x 6 2 (18) 

(dL d dL\ dL 

= Ldt- 1- —-- I (dx - xdt) + —{dx - xdt). 

\ dx dt dx) dx 

Lemma 2.17. For every section cr : / —»1 x R”, and each vector field X = rf + 
on I x R", 

fa*(£ x 30) = fcT*(£ X 2 (Ldt)). 

In particular, if 

f o-*(£ x i(Ldtj) = 0, 

then, 

fa*(£ x 30) = 0. 
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Proof. By proposition [2771 

j 2 cr*£ x i 02 = 0 and j l cr*£ x td l - 0. 

Lifting these equations to / 3 , we get 

/ct*£ X 37tI 2 02 = 0 and j 3 cr*£ x n: 3 ] 6\ = 0. 

Equation ( fT8l) written in terms of pull-backs of the jet projections njj reads 


dL 


where we consider the coefficients - ■j/jj- ') and || as functions on J 3 . Hence, 


j cr*{£ x 30) 


iVfeo;#*))+A* (**= (| - II)) 

+fcr* ^ f(r*(n 3 2 0 2 ) + fcr* fo-*(£ X 3jr* 3a 6 2 ). 


But, 


fcr*(£ x3 (n* 32 (Ldt)) 

fo-*n* 3l 9 l 

j i cr*{n* 32 e 2 ) 

j 3 cr*(£ X 37i* xl di) 

fa*{£ x ,n* 32 6 2 ) 


j 2 cr*(£ X 2 (Ldt)), 

= 0, 
j 2 <T*02 = 0, 

j 3 cr* (77-3 i (£ x i 6»i )) = f ( t *{£ x \ 0 \ ) = 0 , 
f(r*(7T* 32 (£ x 20 2 )) = j 2 cr*(£ X 262 ) = 0. 


Hence, 

fo-*(£ x 30) - fa* (£ x i(Ldt)). 


q.e.d. 

Proposition 2.18. If X = Tj- + fjp is an infinitesimal symmetry of the Lagrange 
form Lclt then, for every section cr of I X R" —> /, 

lV7r = jV(X 3 j0), (19) 

where X 3 /.v the prolongation ofX to J 3 . If cr satisfies the Euler-Lagrange equations 
for L, then Jf<T*(X _J 0) /.v constant. 
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Proof. Omitting pull-backs by j 3 cr for the sake of transparency, we may write 


(X 3 _| 0) = (20) 

d d d d ... ... d 

T — + -1- (t - JCt)— + (f - 2 XT - xf)— + (f - 3 XT - 3xf - XT) — 


dt 


dx 

dL d dL 


dx 


dx 


dx dt dx 


dL d dL 


dL 


_l \Ldt + \ —- —- — (dx - xdt) + —(dx - xdt) 


dx dt dx 


dL d I dL 


dL 


= t\L — x\ —- — ——)-x—\ + f\ — — I + (f-xf)— 


dx 


dx 

dL d dL 


dx dt dx 


dL 


dx 


dL d 


= LT + ^-T,\gjl + = 


dx dt 


Hence, j 3 <r*(X 3 ) _l 0) = j 3 cr* Jfx and is a constant if cr satisfies the Euler-Lagrange 
equations. Suppose that y : I —» / 3 is a section of the source map such that its tan¬ 
gent Ty(I) is contained in ker d®. Since £z@ = Z _l d® + d(Z _l 0), it follows that 
for any infinitesimal symmetry Z of 0, we have 


dy* (Z _l 0) = -y*t/(Z_l 0) = -y*£ z ® = 0, 


which implies that y*(Z _i 0) is constant. This holds for any section y with Ty(I) c 
ker0. Moreover, if y = j 3 cr, where cr satisfies the Euler-Lagrange equations, and 
Z = X 3 , where X is an infinitesimal symmetry of the Lagrange form Ldt, then 
le mm a (12.171) implies that 

y*£ z ® = j 2 a*£ x L dt. (21) 


Hence, if j 2 o-*(£ x ^(Ldt)) = 0, then fcr*(£ x 30) = 0, and j 3 cr*(X 3 _l 0) is constant. 
Moroever, if X = then 


(X 3 _l 0) = (22) 

d d . d d d 

T— + -1- (f - xf)— + - 2 XT - xf)— + (f - 3 XT - 3xf - xf) — 


dt dx 


dx 

dL d dL 


dx 


dx 


= t\L 


dx dt dx. 


dL d dL 
dx dt dx 


t dL d I dL 
,dx dt\ dx 


dL 


_l \Ldt + — —(dx - xdt) + —(dx - xdt) 


dx 


d dL 
dt dx 


(f-TX) + — — (f~TX) 
dx dt 


.dL) 

~ X ~d~x J 

'*«(S 


dL d 


dL 
* dx 


q.e.d. 


Equation ( fl9l ) gives a simple way of finding constants of motion corresponding 
to symmetries of the Lagrange form. 

Example 2.19. If x = (x') are Cartesian coordinates in R", then the action of SO (n) 
on J 3 is generated by vector fields 

d j d d . d ...id d ...,• d ... ; d 
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Hence, for a section cr of 1 x R" —> /, 

U\/ XlJ ) = jVCXj _I 0) = 

= fcr*{x l p x j - x J p x i + x l p xJ - x ] p x i). 

In the following we omit the symbol j 3 cr*, and write 

fxij = x'p x j - x J p x , + x'pp - x ] p X i. 

If L is invariant under the action of SO(3) on J 2 . then jZx, is constant on solutions 
of the Euler-Lagrange equations. 

This suggests that there may be additional conserved quantities coming from 
symmetries of the Cartan form that arc not symmetries of the Lagrange form. 

Definition 2.20. An infinitesimal symmetry of the Cartan form 0 is a vector field 
Z on / 3 such that £z@ = 0. 

Set 

/z = ZJ© 

for each infinitesimal symmetry Z of the Cartan form 0. 

Theorem 2.21. Let Z be an infinitesimal symmetry of the Cartan form and let y : 
I —> J 3 be a section of the source map. IfTy(I) is contained in kcr r/0, then 

yfz = y(Z_1 0) 

is constant. In particular, if cr satisfies the Euler-Lagrange equations, then j 3 a 
is a constant. 

Proof. Since 

£z© = Z_I cl 0 + d(Z _I 0), 

it follows that for any infinitesimal symmetry Z of 0, and any section y : I —* J 3 
of the source map, that 

dy*(Z_\ 0) = —y*(Z _l d@). 

If Ty(I) is contained in kcrr/0. then y*(Z _l r/0) = 0 and y*(Z _l 0) is constant. 

q.e.d. 

2.3.3. Symmetries up to a differential. The Cartan form may yield more conserved 
quantities than those that follow directly from the Lagrangian approach. However, 
even more conserved quantities may arise if the notion of symmetry is relaxed 
somewhat. 

Definition 2.22. A vector field X on I x R is a symmetry up to a differential of the 
Lagrange form Ldt if there exists a function F on J 2 such that 

£ x i{Ldt) = -clF, 

where X 2 is the prolongation of X to J 2 . 
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Proposition 2.23. If a vector field X = r-| + on I x R" satisfies the condition 


£ X 2 (Ldt) = -dF, 

where X 2 is the prolongation ofX to J 2 , and F is a function on J 2 , then 


dL d IdL\ 


dx dt \ dx J 


' x = F + LT + \ — --\ — \Uf-TX) + —-(Z-TX) 


dL d 


dx dt 


is constant along solutions of the Euler-Lagrange equations 

dL _ £(—] + — (—\ - () 


dx dt\dx dt 2 \ dx 


(23) 

(24) 


Proof The pull-back of equation (1231) by the jet bundle projection : J 4 —> J 2 
gives 

n* 4 2 £ X 2{L dt) = -n* A2 dF. 

Therefore, for a section <r of / xR", 


fcr*(nl 2 £ X 2(Ldt)) = -f(T*{n A 2 dF) = -j 2 a*dF. 

Taking into account the form of the Noether identity given in equation (fT6l ) yields 


•2 * r j-i »4 * 

-j cr dF = j cr 


. A * 

+J o- 


dL d ldL\ d 2 tdL 


Tx~d, Ur + dfi te ll(f_Ti)l+ 


dt 


dL d IdL 


dx dt \ dx 


dL d 


~r\L t + — - ~r -rr (^ - tx) + — —{£ - tx) 


dx dt 


If cr satisfies the Euler-Lagrange equations, then 




J cr 


dx dt\dx] dt 2 \dx 


and consequently 


. IdLdldL\\ dL d 

= F + 1 Z.r + - - - - - rt) + --g - ri) 


is constant along cr. 


q.e.d. 


Example 2.24. Probably the most well-known example of this sort of behaviour 
occurs in the first-order theory as the Runge-Lenz vector in the Kepler problem. In 
this case there is a symmetry of the dynamical system that is not lifted from the 
configuration space, and this implies that the Lagrangian is not invariant under the 
action of the symmetry group, but changes by a total derivative. This is discussed 

inlUl. 


In a similar way, infinitesimal symmetries up to a differential of the Cartan form 
are defined as 

Definition 2.25. A vector field Z on ,/ 3 is a symmetry up to a differential of the 
Cartan form 0 if there exists a function F on / 3 such that 


£ z 0 = -dF. 


(25) 
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Proposition 2.26. If£z 0 = -dF, then for a section y \ 1 —> J ' such that T(y(I )) 
is contained in ker <70, the function F + (0, Z) is constant along y. In particular, 
F + (0,Z) is constant along the jet extensions of sections it of I X R" that satisfy 
the Euler-Lagrange equations. 

Proof Since £z© = Z _I d 0 + d (0, Z), equation (1251) gives 

Z_l <70 = -d(F + (0,Z». 

Hence, 

y* (Z_I r/0) = -y*d(F + (0,Z» = -dy*(F + <0,Z». 

However, y* (Z_l <70) = 0 if T(y(I)) is contained in ker <70. Therefore, dy*(F + 
(0, Z)) = 0, which implies that F + (0, Z) is constant along y. 

If a section cr of I x R" satisfies the Euler-Lagrange equations, then T(firil)) is 
in the kernel of <70. Therefore, F + (0, Z) is constant along fir. q.e.d. 


2.4. Parametrization invariance. Let Diff + R be the group of orientation pre¬ 
serving diffeomorphisms of the real line R. Then, for every p e Diff + R, pit ) > 0 
for all t . Each <p e Diff + R gives rise to another diffeomorphism 

(/9°:RxR"-»RxR": ( t, x ) i-» p°it, x) = (<p(t), x). 

The prolongations of p° to jet bundles can be written as follows 


<P 


<P~ 


<P 


0 —> J 1 : ( t, x, x) i p l (t, x, x) = (i pit), x, -], 

\ <Kt)) 


(26) 


2 2 2 / JC JC 

J -» / : (t, X , x, x) p (t, X , x, x) = p(t), X , ——7 - x—- ., 

\ <p(t) ‘Pity) 

J 3 —» 7 3 : (f, x, x, x, x) i-» p 3 (t, x, x, x) = 


, * * . <P 

= \< p ( t ), X , —T7-X 


* ,.. ^(0 ,,( 3(0 , O ,.( 3 ( 0 2 

24 22 , 22 , I .3 22 


0(0 ’ <£(0 2 <£(0 3 ’ s3(0 3 <j3(0 4 ^(O 4 ' s3(0 5 


Proposition 2.27. For p e Diff + R, 

^*01 = 0i, 

p~*6s = —6i. 

<P 


Proof. Let / = (- 00 , 00 ). The contact forms are 0\ = dx - xdt and O 2 - dx - xdt . 
Since 

P l ( t , X , X ) = ( pit ), X , x / p ( t )), 

it follows that 


p u 6\ = dx — x'dt = dx --— dip{t)) = dx - —pit)dt = dx — xdt = 6\. 

ipit) pit) 

Similarly, 


2 . j2 


p~ : J 


J 2 : it, x, x, x) i-> p 2 it, x, x, x) = 


.. X X 
pit),x,—,—— 
pit) pity 
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implies 


<p 2 *0 2 


= ip 2 *(dx' - x" dt) 

■ 

dx x , , , 
= —— ~ ——7^(0 dt 


- x-^—^ | d(<p(t)) 


Pit) 3 
x 


Pit) 

dx 

Pit) 


Pit) 2 
x 


Pit) 2 <p{t) 3 


x———x | <p(t) dt 


X ip(t) 

-?ip(t) dt -— dt + x -- dt 

m 2 w ) 2 


——(dx - xdt) = — d 2 . 
‘Pit) <fit) 


q.e.d. 


A one-parameter subgroup (p E : 1 1 -» t - 1 p E (t) of Diff + R is generated by a vector 
field X T = T(t)d t , where 


T(t) = 


dtp E {t) 


ds 


e=0 


is an arbitrary smooth function on R with f(t) + 0 for all l e R. The Lie algebra of 
the group Diff + R is the collection of vector fields 


diff + R = { X T = r(t)d t \ r e C°°(R), and f{t) + 0 for all t\ 


with the Lie bracket 


[Ti(t)d t ,T 2 (t)d t \ = (ri(f)f 2 (f) - T 2 (t)ti(t))d t . 

Definition 2.28. The variational problem with the Lagrangian L is parametrization 
invariant if the Lagrange form Ldt is invariant under the action of Diff + R on J 2 . 

Remark 2.29. Suppose that the Lagrange form Ldt is Diff + R-invariant. This 
implies that for X T = T{t)d t , the Lagrange form Ldt is invariant under the one- 
parameter subgroup of Diff + R generated by X T . By Theorem 12. 1 81 r fx T — (0, X S T ) 
is constant on solutions of the Euler-Lagrange equations. 

Theorem 2.30. If the Lagrange form Ldt is Diff + R -invariant, then 

jVA = 0 (27) 

for all X T e diff + R and all solutions cr of the Euler-Lagrange equations. 

Proof Recall that 0 = p x dx + p x dx - H dt. Omitting pull-backs by per for the 
sake of cleanliness, 

Jx T = <0,2f 3 > = [pdx + pplx - Hdt, (28) 

= - (Px, x) t - Hr. 

If cr satisfies the Euler-Lagrange equations, then P<t* JLx t is constant. 
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Take two points, to < t\ in I, and consider two other vector fields X Tl and X T2 in 
diff + R such that 

r(t 0 ) = rife) = r 2 fe) and t(t 0 ) = fife) = t 2 (t 0 ), 

T(h) f ri(fi) = r 2 fe) andffe) = fife) ± f 2 fe), 

rfe) f Tife) = r 2 fe) and ffe) = fife) f t 2 fe). 

Then, fefefe, ,/x T] (0 and JPx T2 (f) are constant along / 3 <x. Moreover, the assump¬ 
tion that rfe) = rife) = r 2 fe), Tfe) = fife) = f 2 fe) and equation (1281) imply 
that fxSt) = /x T U) = /xfeO for all t. Therefore, Jx T< (0 - ,/xfe) = 0 and 
^x T ^f) ~ = 0 for all t. Using equation (|28T ) and setting t — t\ yields 

p x fe)xfe)fife) + H(ti)Ti(h) - p x (h)x(ti)t(ti) -H(h)T(t\) = 0 

Px(h)x(h)t 2 (h) + H(ti)T 2 (h) - p x (h)x(h)ti(h) - tffefefe) = 0. 

Since r(tj) ^ rife) and f(q) = fife), the first equation above yields //fe)(rfe) - 
t i fe)) = 0, which implies that //fe) = 0. Similarly, the assumption that rife) = 
r 2 fe) and fife) f f 2 (fi) together with the second equation above yield p x (t\)x(t\ ) = 
0. Since t\ is an arbitrary point in 1 different from to, it follows that 


Hit) = 0 and p x (t)x(t) = 0 for all tel. 
Substituting this result into equation 1281 gives (f27l) . 


(29) 

q.e.d. 


Remark 2.31. Equations (l29l) . rewritten in terms of the configuration variables read 


x 




dL 

X dx 

.dL T 
+ x -— - L 

ox 


0 , 

0 . 


(30) 


These are the identities for our reparametrization invariant Lagrangian that follow 
from the second Noether theorem f lfTTl ). The proof of Theorem |2.30| cstablishcs the 
equivalence between the Noether identities ( l30l) and the vanishing of the constant 
of motion o fx T corresponding to every X T e diff + R. 


2.5. Arclength parametrization. Denote by {x,xf) the Euclidean scalar product 
and by \x\ = xj(x, x), the corresponding norm in R". For a curve c : I —> R ":/h 
x{t), where I = [to, h I, the arclength of the section of c from to to t is 

s(t) = I" life)!^- (31) 

Jto 

In geometric problems it is often convenient to parametrize a curve in terms of its 
arclength. If t is the arclength of c, then along c 


(x, x) = \x \ 2 = 1, (32) 

(x, x) = 0, (33) 

(x, x) + (x, x) = 0, (34) 

<x,x) + 3<x,x) = 0. (35) 
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These equations determine submanifolds M l , M 2 , M 3 and M 4 of J 1 , J 2 , J 3 and J 4 , 
respectively. 

Proposition 2.32. Let X — Tj - be a vector field on the configuration space Q. The 
necessary and sufficient condition for its prolongations X 1 , X 2 , X 3 and X 4 to be 
tangent to M 1 , M 2 , M 3 and M 4 , respectively, is that the restriction ofr to M 1 , M 2 , 
M 3 and M 4 , respectively, is constant. 


Proof. The vector field 


X' = -rf— + 4 


dx' dt 

on J ] , generating the reparametrization transformation, differentiating (x, x) gives 

d d 


td — + T—j (x, X) = -T (x, X) . 

Thus, X 1 is tangent to M l if and only if t = 0. Moreover, 

Next, 

( d . d d 

t— ~ ( 2 *i f + I *> 

dt dx 1 dx' 


and 

Further, 


= -3t(x, x) - f(x, x) 
(X 2 {x,x)) 


M 2 


= ~ T . 


X 3 ((x, x) + (x, x)) = 

( Q Q Q Q 

T — - TX l -: - (2x,t + Xif ) : - (3 Xif + 3x,f + Xif ) : I ((x, x) + (x, x)) 

dt dx' dx' dx' 1 

- -4 f (x, x> - 3r (x, x) - t (x, x) - 4t ((x, x) - 2f (x, x) , 


and 


X 3 «x,x) + <x,x»| M , = 


-T. 


Finally, 

X 4 (<x,x (4) > + 3<x,x>) 

T dt ~ ^~df ~ + jC ‘ f) 'dX ~ + + ^df) ((*’ x(4) ) + 3 + 

- (4x (4) t + 6xr + 4xf + xr (4) ) ^ (4) (^x, x (4) ^ + 3 (x, x)) 

= -5t(x, x (4> ) - 9f(x, x) - 9f(x, x) - 15f (x, x) - 7r (x, x) - r (4) (x, x). 

Therefore, X 4 ({x, x (4) ) + 3(x, x)J = 0 if f = 0 and 


X 4 ((x, x (4) > + 3(x, x>) 


M 4 


= —f 4 \ 


q.e.d. 




ELASTICA AS A DYNAMICAL SYSTEM 


21 


Suppose a local section cr of Q with domain I c R and with j l <r(I ) not in M 1 , 
satisfies x{t) f 0 for all tel. 

Lemma 2.33. There exists p e Diff + R such that 

= 1^1 
at 

for all tel. 

Proof. This follows from the fundamental theorem of the calculus. q.e.d. 


Then, 


and it follows that 


dx dx dt dx 1 x(t) 

dp dt dp dt |x(f)| \x(t)\ ’ 

dx 

dp 


Thus the new parametrization gives rise to a section p*cr with its first jet in M l . 
Similarly, the k-jet of p*cr is in M k . 


2.6. Hamiltonian formulation. The Liouville form on the cotangent bundle T*J ] 
with variables (t, x, x, p ,, p x , p x ) is 

6 = p t dt + p x dx + px dx. (36) 


The exterior derivative 


u> = d6 


is the canonical symplectic form of T*J ] . 
Lemma 2.34. The action 


Diff+ Rxi 1 

lifts to an action 


J 1 


{p, (t, x, x)) 



(37) 


Diff+R xl'/ 1 -» T*J l : (38) 

(p, (t, x, x, p t , p x , pf)) i-» (p{t), x, p(ty l x, p t p(ty l + (p x , x) p(ty 2 p(t), p x , p{t)p x ) ■ 

The lifted action IvS’D is Hamiltonian with momentum map JT : diff + R* 

such that, for X = r{t)d t £ diff + R, 


x, x, p t , p x , p^ = T(t)p, - fit) (p x , x). (39) 

Proof. The action of Diff + Ron/ 1 is given by 

Diff + R x J 1 —» 7 1 : ip, it, x, .i)) i-» it', x', x!) - [pit), x, pity l x}, 

(see equation (l26l).) The lifted action takes ( p, , p x , p x ) to (/;', pf pf such that 
p'jdt' + p' x dx + p\dx = p t dt + p x dx + p x dx. 
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But dt' = ip(t)dt, dx' = dx and dx' = ip(t) 1 dx - ip(t) 2 xip(t)dt. Hence, 

p' t dt' + p' x dx' + p' k dx' — p' t Cp(f)dt + p x dx + p' k <p(ty x dx - {p\, x) ip(t)dt 
= (p',<p(t) - (p k , x) <p(t)~ 2 ip(t )) dt + p x dx + p' k (p(t)~ x dx 
= p t dt + p x dx + p x dx. 

Hence, 

p'M* ) - (p' x , x) <p(ty 2 (p(t) = p t , 

Px' = Px, 

P k ip{ty x = p x . 

Therefore, 

Pt = Pt<p(t)~ X +(Px,x)‘p(ty 2 ip(t), 

Px = Px, 

Px = <P(t)pt. 

The lifted action on the cotangent bundle is Hamiltonian, and the value of the 
momentum map on an element of the Lie algebra is given by the evaluation of the 
Liouville form on the vector field generating the action of the one-parameter group 
corresponding to this element of the Lie algebra. Hence, for X = T(t)d, e diff + R, 

xf X {t, X, X, Pt, p x , Px) = , X) (t, X, X, p t , p x , p x ) = (e,X xS }(t,x,x,p t ,p x ,p x ) 

= [p,dt + p x dx + Pxdx, r-^ - 

= T(t)p, - f(t) {p x , x) . 


q.e.d. 


Definition 2.35. The Legendre-Ostrogradski transformation 

2z? : 7 3 —* T*J X : (t, x, x, x, x) i-» (t,x,x,p t ,p x ,p x ), (40) 

is given by 

Pt = ~H = -xp x - xp x + L 
dL 

Px = -777; 

ox 

dL d dL dL d 

dx dt dx dx dt^ x 

The Legendre transformation was extended by Ostragradski in lfT8ll . However, 
for brevity, from now on we shall just refer to it as the Legendre transformation. 

Clearly, 2z? is a smooth map of into T*J X . If .if is a diffeomorphism, Ostro- 
gradski’s approach leads to a regular time-dependent Hamiltonian theory with the 
Hamiltonian 


H = -p t = xp x + XPx - L. 
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In geometric problems, the Lagrangian is often reparametrization invariant. In 
this case L does not depend on t and the range of the Legendre transformation is 
restricted by the equations (see (12.301) 1 

H = 0 and (p x , x) = 0. 

Theorem 2.36. If the Lagrange form Ldt is invariant under the action of Dili’, R, 
then the Legendre transformation S£ intertwines the actions of Diff + R on J 3 and 
on T*J l . 


Proof The map «£?, defined by (l40l) . intertwines with ip 3 , the induced action of 
Diff + Ron/ 3 , if and only if 

ip 1 o J£(t, x, x, x, x) = 2z? o tp 3 (t, x, x, x, x), (41) 


for every (t, x, x, x, x) e J 3 . Here the map <p l : T*J l —> T*.J l is the lifted action of 
<p l to the cotangent bundle T*J l , defined by (1381) . 

Let us first compute the right hand side of (|4TT) . Let L': J 2 —> R be a function on 
J 2 such that ip 2 )*L' = L. This function is an induced Lagrangian by the Diff + R- 
action. Based on (l40l) . the map «£? o tp 3 can be calculated by: 


it, x, x, x, x) ^ (t', x, x',p',p' x , p k ) 
t' = (pit) 

, 4 

x = — 

<P 

Pt = -(Px, x ')-(p'x’ X, ') + L ' 
, _ dl! 

P *-~d7' 

, dL' d , 

= 7 - izW 


x = 


dx' dt' 
x ipx 

<p- (pf 


In order to calculate the partial derivatives of I!, we use the assumption that the 
Lagrange form Ldt is invariant under the Diff + R-action, i.e., 

(tp 2 )*(L'dt') - (L' o (p 2 )ipdt = Ldt => L = Cpil! o p 2 ). 


Therefore, the partial derivatives of L' with respect to x, x' and x" are calculated as 
dL _ . dL' dV _ 1 dL 

dx ^ dx dx tip dx 

dL _ . IdL' dx" \ _ 1 dV dl! _ . dL 

dx ^ \ dx" dx) ip dx" dx” ^ dx 

dL . / dL' dx' dL' dx" \ dL' ip dL' dL' dL ip dL 

dx ^ \ dx' dx dx" dx ) dx' ip 2 dx" dx' dx ip dx 


Since Ldt is parametrization invariant, by theorem (12.301) . L is time invariant, p, = 
-H = 0 and (p k , x) = 0. As the result, we show that the partial derivative of L' 
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with respect to t' is equal to zero: 

„ dL T . IdL'df dL'dx' dL' dx" 
a * \df dt dx' dt dx" dt 


dt 

- p(l' + — 

\ P 


dL' \ 2 1 dL’ \ 3 ip I dL' 

-, jc I-I - V ) 4-( - T 


dx' 


P 


dL 


dx 

dL 


P 


ip \dx r 


,x) + 


ip 2 \dx‘ 

3 <P 


;,x 


. 9 dL' ip I dL' 


= “| L \ (P*’ jC ' > ~ 2 ( pi ' ^ + ~2 (P*’ I + - - (P*’ ^ 


p 


p- 


+ tfi ~dr p 2 \dx" ,x 
. 9 dL' ip 


P 




p 


dx 

dL 


dt' 


= -\ L -\-qI- - f {Px^) + (px,x)- 2 (p,,x)) + ip — 


dt 


■ 2 


dL' 


dt' 


ip n , .. , .... .idL' ip , 2 dL' . 2 dL' dL' 

= -{L- < p x ,x > - (p it x)) + p-— = -p, + p — = pr— => -p— = 0. 

p dt' ip dt' dt' dt' 

Now, we can calculate the right hand side of (HTl) in terms of elements of 7 3 : 

, dL . dL 

dx 

dL 

dx 


Px = 


= P ^77 = PPx 
dx 


dL 


, ^ d , dL ip dL Id i.dL\ ^ 

Px = ~^~ dF iPi) = dl + fdX ~ pJt \~dic) = dl~ Jt {Pi) = Px 


p' t = -(p' x ,x')-(p'„x") + L' = 


x 

Px, ~ 
P 


X px 
PPx, -T - -T ) + T 


r p j 


p 


= -7 (- (Px, X) - (Px, X) +L) = -Pt - 0, 
P P 


where p x and p x are considered as functions on J 3 , defined by: 

dL d ldL 


Px dx dt \ dx J ’ 
dL 

Px = "7777 ■ 
dx 

In the following, we calculate the left hand side of (I4TT) . Since Ldt is Diff + R- 
invariant, and based on (|38T ): 

p l ° ■Z'it, X, X, X, X) - p\t,X,X,0,Px,Px) = (p(t), x, x/ip, 0, p x , PPx)- 

Therefore, for Diff + R-invariant Lagrange forms the relation (|4TT) holds. This 
completes the proof of the theorem. q.e.d. 


Corollary 2.37. If the Lagrange form Ldt is Diff + R -invariant, then the Cartan 
form 0 is Dilf + R -invariant. 

Proof Since Ldt is Diff+ R-invariant, theorem (12.36b implies that for p 6 Diff + R, 
p 3 *££* - L£*p u , where p x is the lift of p ] to the cotangent bundle T*J { . But, 
0 = 2z f*6, and the Liouville form 6 is invariant under the the lifted action p L 
Therefore, 


p 3 *@ = p 3 *sro = jgf *p u e = sro = 0. 
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q.e.d. 

Corollary 2.38. The range of the Legendre transformation S£ : ./ ' —> T*.f is 
contained in the zero set of the momentum map J? : T*J l —> diff + R*. That is, 

J£(J 3 ) C 


3. Classical elastica 


3.1. The variational equations. The elastica functional is given by 


A[cr] = 


/'V 

Jtc, 


\x\ dt, 


where 

cr : [t 0 , t\\ —> [to, fi] x R’ : / h (t, x(f)) 

corresponds to a curve cr : r i—> x(t) in R 3 , and k is the curvature of cr. Since 
the curvature of the curve depends on its second derivatives, this is naturally a 
second order variational problem. As the curvature in Cartesian coordinates x = 

(x 1 ,* 2 ,;*: 3 ) e R 3 is 

|x| 2 (x, x) 


k 2 = 


l *| 4 


1*1 


• 16 


the elastica Lagrangian is 


L(x, x, x) = 

1*1 


<*, x ) 2 

l *| 5 


(42) 


It is defined on {(t, x, x, x) e J 2 \ x + 0). 


Proposition 3.1. The elastica Lagrangian 442 D is invariant under translations and 
rotations in R ’ and is independent of parametrization. 

Proof The expression (I42T ) for L is independent of x and depends only on Eu¬ 
clidean scalar products of x and x. Hence, L is invariant under translations and ro¬ 
tations. Moreover, the curvature k of a curve is independent of its parametrization, 
and |x| dt = ds is the element of arclength. Therefore, Ldt = k 2 ds is independent 
of parametrization. q.e.d. 


For elastica, Ostrogradski’s momenta are 


and 


Px = — 


(x, x > 5/2 

The Euler-Lagrange equations 


„ X „ (x, x) X 
Px = 2— - 2AAJ_, 
|x| 3 |xf 


., ... ,. .... < x,x)x (x,x>x 

«X, x) X - (x, x) x)-+ 0- 


(x, x) 


l\5/2 


(x, x> 


v\5/2 


c (x, x ) 2 X 
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can be written in the form 


dL d 

& - j , p ‘ = °- 


Since the Lagrangian is parameter-independent, the Euler-Lagrange equations are 
necessarily degenerate in that they do not determine the fourth derivative x uniquely. 
Let 


jc*l = ^$x and x ± - x - $ 

|i| 2 

denote the components of x that are parallel and perpendicular to x, respectively. 
The Euler-Lagrange equations written in terms of this decomposition are 


jc~ L 


= 6 \x\ 2 (x, x) 3c+4 (x, x) x + - |x| z x - 10 (x, x) (x, x) - - 


5 |x| 2 (x, x) . 


I-^I 


■ i2 


35 (x, x) 

x+— - —x 


4 lil e 


(43) 


They determine x x , but leave the component undetermined. On the other hand, 
the parametrization-invariance of the problem allows us to use parametrization by 
the arclength. In the following, assume that t is the arclength parameter of the 
curve. Therefore 


\x\ 2 - (x,x) = 1, 

(44) 

and, by differentiation 

(x, x) = 0, 

(45) 

{x, x) + (x, x) = 0, 

(46) 

(x, x) + 3 (x, x) = 0, 

(47) 

as well. Substitution into (1431) and (1471) yields 

3 

x ± = -- x| 2 x and - -3 (x, x) x. 

(48) 


These equations determine the elastica completely. In other words, the choice of a 
parametrization determines an equation for the component 


Remark 3.2. This yields an equation of the form 

x = fix, x, x, x) 

to which theorems in differential equations apply that guarantee the local existence 
and uniqueness of solutions. 

3.2. The Frenet frame. The elastica equations (1481 ) are conveniently studied in 
the moving frame ( T , N, B), where T = x is the unit tangent vector, N the normal 


vector and B the binormal vector of the 

curve 1 1 —> x(t). The Frenet equations 

are 

T = 

kN, 

(49) 

N = 

-kT +tB, 

(50) 

B = 

-tN, 

(51) 


with k = \x\ the curvature and r the torsion of the curve. In order to relate the 
torsion r to the derivative variables, observe that 

x-'kN + kN-kN- k 2 x + ktB, 


(52) 
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which implies that, if k ± 0, 

r = k~ 1 (B, X > = K~ l (T X N, X) = K 1 (x X X, x>. (53) 

Differentiating (1521) and the Frenet equations imply 

x — -3kkT + (if - k 3 - kt 2 )N + (2a :t + kt)B. (54) 

This, together with equation (|48] | implies that 


2'kt + kt - 0, (55) 

2k + k 3 - 2kt 2 = 0. (56) 

Equation (l48l ) does not lead to any new condition because k = |x| implies that 
(x, x) = k'k. Equation (l55l) can be immediately integrated to yield 

k 2 t = c, c a constant. (57) 


If k ± 0, substituting r = 4 into equation (1561) and integrating gives 


\ ^2 

k 2 + -k 4 + — = constant. (58) 

4 k- 

Integration of equation (l58l) determines completely the functions k(i) and r(t) in 
terms of the initial data k(?o)> ^(Xi) and r(to)- Thus, in order to find the solution 
t h-> x(t), it suffices to integrate Frenet’s equations assuming that the curvature k and 
the torsion r are known functions of t. This can be achieved using the conservation 
laws for elastica. 


3.3. Conserved momenta. Since the Lagrange form for elastica is invariant under 
translations, it follows that the linear momentum p = jr - Jj {'jp. ) is conserved. In 
the arclength parametrization 

- 2(x - (x, x) x) - (x, x) x = p = constant. (59) 

The arclength parametrization implies 

p = -2x - 3 (x, x) x. (60) 

Similarly, rotational invariance of the Lagrange form implies that the angular mo¬ 
mentum 

Aj = X l p x i - x J p x , + x‘p x j - X J p x , 

is conserved (see example [2.191 ) Setting I - (l\,h, h), where 

U = tijk^Xjk, 

gives 

l = x x p x + x x px- (61) 

The expression (1421) for the elastica Lagrangian in arbitrary parametrization yields, 
in the arclength parametrization 

Px = 2x - x, 

which implies that 

1 = x x p + 2x x x. (62) 
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Proposition 3.3. The conserved momenta in the moving frame ( T , N, B ) are 

p = -k~T - 2kN - 2 ktB , (63) 

I = x x p + 2kB. (64) 

Proof. Equation (1521) implies that 

p = -2(kN - k 2 T + ktB ) - 3k 2 E = -/rJ - 2/dV - 2ktB. 

while (l49l) and (BTT) yield 

xxx=Tx ( kN ) - /dT x N) = kB. (65) 


q.e.d. 


Proposition 3.4. Scalar equations for the curvature and torsion ( 1571 ) and ( 1581 ) can 
be rewritten in the form 


, 1 

kt = --(Ip), 

■ 2 1 4 1 , ,2 

* + 4* + l6^ = 4 IPl ' 


Proof Equations (1631) and (l64l) imply 

(l,p) = 2 k(B,p) = 2,k(b, -icT - 2kN - 2 ktB^ = -4a: 2 r. 
Equation (l63l) implies 


,AJ2 


\pf = rf + 4A: 2 + 4a: 2 t 2 - 4(k z + -k 4 + 


•2 


= 41 k 2 + 4/c 4 + 


1 


(IpY 

\6k 2 ] 


( 66 ) 

(67) 


q.e.d. 


Suppose that a: and r as known as functions of t (the differential equations im¬ 
ply that they may be expressed as elliptic functions.) It remains to show how the 
integration of the Frenet equations is aided by the conservation laws 

p = -2kN - 2 ktB - k 2 T, 

1 = xx j r + 2kB. 


Theorem 3.5. The velocity of the elastica in the direction of the consented momen¬ 
tum p is 


(X,p) = -K 2 . 


( 68 ) 


Hence 


(x, p) = (x(to), p) - 



KTds. 


(69) 
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Proof. Taking the scalar product of / and p yields 

< B,p > = 

2k 

Moreover, by taking the derivative of both sides of this equation and using the third 
Frenet equation, 


7, (B. p) = (B. p) = -r (N, p) = | (T) «, p) = XL <(, p >. 


Therefore if r ^ 0, by (1661) we have 


Finally, the second Frenet equation yields 


</,p> = -2k. 


— (iV,p> = -K(x,p) + t(B,p) = -2k. 
at 


Therefore if /e + 0, (l56l) implies 

, , 2k T , -K 3 + 2 KT 2 

(x, p) = — + - (B , p) = - 

K K K 


+ ■—- (l, p) = -k 2 + 2t~ + —^(-4 k 2 t) = —k. 


2k 1 


2k 2 


In particular, 

(x,p) = (x(t 0 ),p) + 


f (Ms), 

to 


p) ds = (x{to),p) 


f 

to 


K 2 ds. 


q.e.d. 


It remains to determine the component of the motion perpendicular to p. If x is 
not parallel to p, then the vectors xx p and (x x p) x p span the plane of directions 
perpendicular to p. Their lengths are 


and 


\x x p\ = \-2kB + 2 ktN\ = xj4k 2 + 4 k 2 t 2 = -\j\p\ 2 - k 4 , 
\(x x p) x p\ = ^\p\ 4 - | p\ 2 if = \p\ xj\p\ 2 -k 4 . 


Let D and E denote the unit vectors in the direction of xx p and {xx p) x p, 
respectively. Equations (1631) and (1641) give 


D - 


xx p 

\x X p\ 


= (lp| 2 -^ 4 ) 1/2 (~2kB + 2 ktN) . 


(70) 


Similarly, 


E = 


{xx p)x p —{4k 2 + 4 k 2 t 2 )T + 2k 2 kN + 2k 3 tB 
\{x x p) x p\ 


\p\{\p\ 2 -K 4 ) 


1/2 


(71) 
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Proposition 3.6. The frame (JXt), E(t)) satisfies the equations 

b = - {l ’ Pm E, 


E = 


2(\p\ 2 -K*) 

{Up) \p\ 


2(lpl 2 - K 4 ) 


D, 


where k 2 t = -\{l, p)- 

Proof. Equations (1631) and (l64l) give 

x x p = —2kB — 2kB = —2kB + 2 ktN 

and 

(x x p) x p = (p, p) x - (x, p) p = (4A- 2 + 4/Tr 2 )r - 2 a- 2 kN - 2k 3 tB. 
Differentiation yields 
d 

dP xp> 


(72) 


and 

d r 

dP 


XX p 

= kN x(- 

-2k 2 tT + k 3 B 


2 K-T 

(4k 2 

+ 4 a- 2 t 2 ) 


2k- t 

(4k 2 

+ 4 a- 2 t 2 ) 


2 K-T 

(4k 2 

+ 4k 2 t 2 ) 


2 K-T 

(4k 2 

+ 4 a- 2 T 2 ) 

= (x 

X p)x P 


{(.i x p)x p + 2k 2 kN + 2a- 3 tB) + ic’B 


{x x p) x p 
{xx p)x p 


4a- 4 ATT 


-N + 


4k 2 


(4A 2 + 4a- 2 t 2 ) (4a: 2 + 4a: 2 t 2 ) 

2 k 3 k 

(2 ktN - 2kB) 


k 3 B 


(4k 2 + 4 k 2 t 2 ) 

2/c 3 ^ 

{ixp}xl ’-WTJtvP xph 


2 k 2 t 2k 3 k 

(xx p)x p - —-— (x x p)\x p 


(4k 2 + 4a- 2 t 2 ) 
2 k 2 t 


(4k 2 + 4a: 2 t 2 ) 
2 k 2 t 


((x X p) X p) X P - 


(4 k 2 + 4a: 2 r 2 ) 
2 k 3 k 


(4A 2 + 4a: 2 t 2 ) 

.3 


(xx p) X p 


- ( tf + 4,W ~ m2 + p {i • p))xp ~ xp)xp 

2k 2 t\ P \ 2 \ 2k 3 k 

xx p ——-— ^ (x x p)x p. 


(4k 2 + 4a- 2 t 2 ) (4k 2 + 4 k 2 t 2 ) 

Now compute for the orthonormal frame 

j|i x p\~ ] xx p, \(x x p) x p\~ x (x x p) x pj. 


Since 


1-4 


x X p\ — | — 2a :b + 2at«| = V4A: 2 + 4k 2 t 2 — ^p'p - a 4 , 
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it follows that 


— \x x p\- 1 = —(|p| 2 - k 4 )“ 1/2 = -~(\p\ 2 - K 4 r 3 , 2 (- 4 K 3 k) = 2 k 3 k{\p \ 2 - K 4 )- 3 ' 2 . 
at at 2 

This further implies 


d ,. , 2 k-t . 

7, (XXp> = ~ (^ + 4,W XXP)XP 


2 k 3 'k 


\p\- - A 4 
d ll?T |p 2 | 

S Uxp)x " = W7777) 


( 4 k 2 + 4 k 2 t 2 ) 

2 K-T 2 k 3 k 

- - -(x xp)xp- —- -(x x p), 


(XX p) (73) 


XX p 


2k 2 t\ P 2 \ 

— 2 - -7* XP~ 7~a - 7 (* xp)x p. 


\p\ - K 4 
2 k 3 k 

( 4 k 2 + 4 k 2 t 2 ) 
2 k 3 k 


(x x p) x p 


\p\ - K 4 


\p\ - K 4 


Similarly, 


|(ixp)x/3| 2 = \-x\p\ 2 + p(x,p)\~ -\p\ 4 -\p\ 2 K 4 , 


and thus 


^-\(xxp)xp\- 1 = ~(\p\ 4 — lp| 2 ir 4 )- 1/2 — |p| _ i ( I/?I~ — a'4 ) - 1/2 

clt at at 


2k 3 k 

\p\ 


(\p\ 2 - K 4 r 3 ' 2 . 


Therefore, 


— (\xxp\ 1 * x p) = —(\xxp\ 1 )x x p + \x x p\ l —(xxp) (74) 
at at at 


2 K-T 


(ipl 2 - k 4 ) 

2 k 2 t\p\ 


3/2 


(xx p)x p 


(|p| 2 - /r 4 ) \(x xp)xp\ 


(xx p) x p. 


(75) 

(76) 


Similarly, 

^(\(x x p) x p| _1 (x x p) x p) = |(x x p)x p| _1 j (xx p)x p + 


+ |(d x p) x p | —((x xp)xp) 
at 

= (|p| 2 - k 4 )~ 3 / 2 2 k 2 t \p\ xx p 
2 k 2 t\p\ \x x p\ . 


(|p| 2 - a 4 ) 3 / 2 I* x p| 
2 k 2 t\p\ 1 

(|p| 2 - k 4 ) \x x p\ 


XX p 


xx p. 
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Thus, 


D = 


E = 


(ipl 2 - k 4 ) 

2/c 2 t |p| 


(lp| 2 - K 4 ) 


D, 


and the proof is finished since /rr = (/, p). 

Define the curve of complex-valued vectors Z(t) by 

Z(0 = D(t) + iE(t). 

Proposition 13.61 implies 

Z = D + lE = -! -r--Z. 


Proposition 3.7. Define 

= (Up)\p\ 

and set Zq = Z(fo) = Zo = Do + /2?o> ffew 
In particular. 


2(lp| 2 - k 4 ) 

f'—l 

Jt n (Ipl 2 - 


(IPI - K 4 ) 


<7.q 


Proof. This follows immediately upon differentiating 


<7 z <7 
<7r dt 

and Z(to) = e“^ ao) Z 0 = Z 0 . 

Note that 


Z = —Z = —e~ m Z Q = -ife~ im Z 0 = -i 


(x x p) x p = - |p| 2 ,i + (x, p) p 


(d p) Ipl 
2(|p| 2 - x 4 ) 


Z, 


implies that 


:= - Ipl ‘(fxp)xp 


is the component of i perpendicular to p. 


Theorem 3.8. The time evolution of Xp is 

(|p| 2 - Kfff' 2 


X D (t) = 


\P\ 


(- sin < p(t)Do + cos (p(t)Eo). 


q.e.d. 

(77) 



N 

II 

1 

S; 

(78) 

D(t ) 

= cos f{t)D (> + sin (f>(t)Eo, 

(79) 

£(0 

= - sin f(t)Do + cos <p(t)Eo. 

(80) 


q.e.d. 


Hence, the component t i—> xo + Xp(t) of the motion of the elastica in the plane 
perpendicular to p through xo = x(tf) is 

7 j r‘ 

xj(t ) = - |p| -2 (xo x p) x p - — (|p| 2 - 4 s) 4 ) 1/2 (- sin <p(s)D 0 + cos <f>(s)E 0 ) ds. 

' \P\ Jtn 
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Proof. 


\(x xp)x p | 

-5- £ 

\p\ 2 

\p\ (\P\ 2 ~ ^ 4 ) 1/2 ^ 

1/2 

—(- sin (p(t)D() + cos <p(t)Eo). 


\P\ 2 

(\p\ 2 ~ * 4 ) 

\pT 


q.e.d. 


Corollary 3.9. The elastica equations in the arclength parametrization 

cr : I —> R n : t i-» x(f), 

/zave a unique solution 

v 1/2 


f 

Jtf) 


x{t) = Ao + 
for initial data in 


k(s) 2 (l p \ 2 - *<r) 4 ) 


-(- sin f(s)Do + cos f(s)Eo) 


\P\ 


Mq = {(/, A, A, X, X) € M 3 | /C ^ 0, 0). 




It remains to consider the special cases when k ± 0 and r = 0, and when k = 0. 
Equation (l66l ). k 2 t = shows that (/,p) = 0. Hence, if either k or r 

vanishes at some point to, then it vanishes for all t for which the solution exists. 

(1) If r = 0 and k + 0, the Frenet equations are T = kN, N = -kT , B = 0, and 
the conservation of the linear momentum p and the angular momentum 1 
are 

p = -2 kN-k 2 T, (81) 

1 = xx p + 2 kB. 

Thus, there is an additional conserved quantity, 

1 

B = —xxx. 

K 

As the scalar product (x, p) = - k 2 (x, x) = - k 2 , 

(x(t), p) = (x(t 0 ), p) - f k 2 (s) ds. 

Jt 0 

Taking the cross product of equation (I8TT) with B yields 
-2 kN x B - k 2 T x B = p x B. 

Since N x B = T and T x B = -N, 

-2 kT + i^N = pxB. 


Therefore, 


(x, p x B) = -2k, 
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and 


pxB) 


(x(to ), pxB) + 


(x(to), p x B) + I (x(s), p x B)ds 


to 


(x(t 0 ), pxB) - 2k(j) + 2 K(to). 


Thus, if p x B ± 0, then 



(2) The special case p x B = 0, p ± 0. If p x B = 0, and p + 0, then p is parallel 
to B , and equation (ISTT) implies that k — 0 , so the solution is a straight line. 

(3) If p = 0, then equation (IbTT) implies that k = 0. If k = 0, then x is constant, 
and the motion is again a straight line. 

3.4. Closed elastica. As mentioned in the introduction, a significant motivatation 
for this work was understanding how symmetry and conservation laws could be 
systematically exploited to integrate the elastica equations. In the case of closed 
elastica, as studied by Langer and Singer lfl2ll . it is necessary to add an arclength 
constraint to the variational problem. This results in studying a modified problem 
with an undetermined Lagrange multiplier. The modifications to our analysis are 
straightforward insofar as the use of the conserved quantities is concerned. How¬ 
ever, since there is also an immediate integration and reduction of order in the prob¬ 
lem, which results in a loss of manifest Euclidean invariance, it seemed preferable 
to avoid the arclength constrained problem and keep the full symmetry in order to 
see more clearly how the conservation laws enabled the integration, which is the 
route taken in the previous section. 

In more detail, a slicker, but less transparent approach to the Euler-Lagrange 
equations runs as follows. For a second order Lagrangian L(x, x, x, t) the Euler- 
Lagrange equations are 



it follows that if 



which is the case in the elastica problem. 



immediately integrates to 



where c is a constant. 

Now put q = x, q = x, and set 


I(q, q,t) := L(x, x, x, t) - c ■ x 
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Then the integrated equations are now the Euler-Lagrange equations for the first 
order Lagrangian / 


dl_ 

dq 


d_ldl\ 

dt\dq) 


= 0 . 


3.4.1. Reduction for elastica. 

Lineal - momentum. The elastica functional for fixed arclength is 

k 2 + A ds. 

Here A is a constant whose value is a priori unknown. The curvature k is 

\x x x\ 



K = 


w 3 


so since ds = \x\ dt , the reduced Lagrangian is 

7/ * \ \ C l X d\ 2 , „ , 

Kq,q ) = . ,c + m\ -c-q 

\q? 

It remains to compute the Euler-Lagrange equations and look at them in the Lrenet 
frame. The derivatives are 




\q X q\ A \ 2 


91 2 7. , 

= —rq X {q X q). 
dq |c/| 5 

Define new vai'iables T, N, B and v by setting v = \q\,T = v^ l q, 
{qxq)x q \q \ 2 q - (q, q) q 


N = 


I q X q\ \q\ 


I q X q\ \q\ 


B = T x N. 


This implies 


, dl 2 

q = vT + v~kN, — = -kN. 

dq v 

If we recall the Lrenet equations then it follows that 

dl o KV 

— = (-3k 1 + A)T - 2—N - c, 


and 


dq 


?L[?L\ = 2\-k a T + \-\ N + ktB 


dt \dq 


This implies that the Euler-Lagrange equations are 


(A - k 2 )T - 2-N - 2ktB = c. 
v 

Taking the inner product of this equation with itself and choosing the arclength 
parametrization (so v = 1) yields 

4 (k ) 2 + (A - k 2 ) 2 + 4 ter 2 - c 2 . 



36 


L. BATES, R. CHHABRA AND J. SNIATYCKI 


Angular momentum. The reduced Lagrangian / is not invariant under the rotation 
group SO(3), but it is invariant under the SO(2) subgroup generated by the vector 


field 


d 

X = (cx q) — 
oq 

which is rotation about the axis defined by c + 0. 
mentum is 


The associated conserved mo- 


j = ( pdq,X) = (2v 1 kN,c X q). 


The only nonzero component of this contraction is in the N component of c x q, 
and since the Frenet frame is orthonormal, 


j = 2 k(c,B). 

Taking the inner product of c with the Euler-Lagrange equations gives 

(c, B) = -2 KT, 

and this implies that the conserved angular momentum j is 

j = -4 k 2 t. 

Thus 4 k 2 t 2 = j 2 /4k 2 , and substituting back into the equation for k' yields 

4 Wf + (A - k 2 ) 2 + = c 2 . 

This recovers equations (3) and (4) of foltinek 0, together with the interpretation 
of c as linear - momentum and j as angular momentum. 


4. ElASTKA AS A CONSTRArNED HAMUJONrAN SYSTEM 

4. 1. Range of the Legendre transformation. In this section we discuss the range 
of the Legendre transformation 

2 ¥ : Jq-> T*Jq : (t,x,x,x,x) ^ (; t,x,x,p t ,p x ,Px) 

for elastica with Lagrangian 

(jr 1 -, x.^) 

L(x, X, x) = -r — . (82) 

| 4| 3 

Here, the subscript 0 denotes that 4^0, and 



(')L A'- 1 

(83) 

Px = 

a* l*F’ 

dL d dL 2 , 6 . (x ■ L ,A ± ) 

„. . = ,.v + <x, X) X X, 

dx dt dx |x| 3 |x| 5 |x| 5 

(84) 

Px = 

Pi = 

L - (p x , x) - (px, x) = 0. 

(85) 


The Diff + R-invariance of L dt is responsible for the vanishing of p, above, and it 
implies that the variable px satisfies the equation 

(Pi, 4> = 0 


( 86 ) 
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(see remark [2.311 ) Note that equations (l83l) and ([84b imply 

, \ 4 r 

Wl> = 


(87) 


(Px,x) = (2— -,x)-2 


l*l 3 ’ 

2 _ 

uf 


(x^,x) Ax^,X^) 


\x\ 


m3 


= 2 


1 


, o =2 L= f-\x\ 3 (p x ,p x ), (88) 

1*1 2 


_ 6 _ 

\M 5 


(Px, x) = (-—x 1 + 7777 (x, x) x 1 - - '"//l ' x, x) = -4 1*1 3 (Px, Px) -(89) 


(.V-..V ) . 


1*1 


1 

4 


Equation (I85T) . written in terms of the variables on T*j\, reads 


H(x,x,p x ,px) = -- |*| 3 (Px,Px) + ^ l*| 3 < Px,Px) ~ ^ l*| 3 (Px,Px) = 0. 

Hence, it does not introduce further restrictions of the variables ( t , a, x, p t , p x , 
Equations ([88b and (f89b lead to the new constraint equation 


(Px, x) + 2 (p x , x) = 0, (90) 

while equation (f90b in terms of the variables on T*J ] {] is 

1*| 3 (p.x, Px) + 4 (p x , *> = 0. (91) 

Theorem 4.1. The range of the Legendre transformation — » T*J^ is the 

common zero set of the three functions p„ (p x , x) and (p x , pf) + (p x , x). That is, 

range «Sf = {(?, x, x, p t , p Xt p x ) eT*J l Q \ p, = (p x , x) = |x | 3 (p x , p x ) + 4 {p x , x) = 0}. 

(92) 


Proof Equations (l85l) . ([86b and ([9Tb imply that 

range «£? c {(t,x,x,p t ,p x ,p x ) £ \ p, = (p x ,x) = |x | 3 ( p x ,p x > + 4 (p x ,x) = 0}. 

Suppose (t , x, x, p t , p x ,Px) £ T*Jq is such that p t = 0, (p x , x) = 0 and |x| 3 (p x , p x ) + 
4 (p x , x) = 0. Since (p x , x) = 0, it follows that p x = pf, and equation (183b implies 
x 1 - = 5 |x| 3 Px- By definition, the vanishing of |x| 3 (p x , p x ) + 4 ( p x , x) is equivalent 
to equation (190b : that is (p x ,x) + 2 (p x ,x) = 0. Splitting equation ([84b into its 
components perpendicular and parallel to x gives 


Px = 


Px = 


2 ..., 6 , 

-T* +— ~{x,x)x , 

l *| 3 l *| 5 

(x^fx^) . 

-x. 


1*1 


m5 


(93) 


Equation (1931) gives 
2 


i*i 3 


x L - 


-Px + ~~Z <*, X) X 1 - = -pf + <X, X> ^ |x | 3 Px 

1*1 1*1 2 




- p f + —2 <*, X) p x . 


( 94 ) 
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where (x, x) is arbitrary. Equation ( [94l) is equivalent to equation (l90l) because 

1 

"2 


/ -X /|| -\ {x »* > . x (x- L ,X ± ) 1 . 

(p x , X ) = (p\, X) = - — r (X, X) = -— = -- (Px, X) 


by equation 1 

The above argument shows that the fibre of f£ over the point (t, x, x, p t , p x ,p x ) € 


n * J 1 
J ( 
fact, 


T*Jq such that p t = 0, (p x , x) = 0 and |d| 3 (p x , pf) + 4 (p x , x) = 0 is not empty. In 


where 


2Z? l (t, X, X, Pt, Px,Px ) = 1 it, X, X, x ± + x\x ± + jc*'}, 
X 1 - = - |x| 3 Px, 


q.e.d. 

t 
o- 


x 1 - = 1 (- W 3 Px + 3 \x\ (x, *") px), 

and x " and x^ are arbitrary. 

Theorem 4.2. The range of the Legendre transformation is a submanifold of T* . 
Proof Fix (t, x, x) e J x {] . The constraint equations 

(Px, x) = 0, 

{x,x) 312 (px,Px) + 4(p x ,x) = 0, 

give 

p\ = 0, 

pi = -^(x,x) 3/2 (Px,Px) = ~^(x,x) 3/2 (pf,pf) . 

By assumption, x + 0, which implies that the splitting of vectors into components 
parallel and perpendicular to x is smooth. Therefore, 

{(f, x, x, p t , Px,Px) e T*7 q I p t = 0, (p_x, x) = 0 and | x| 3 (p x , pf) + 4 (p x , x) = 0} 
is equal to 

!( t,x,x,p t ,Px.Px) e T*Jq I p, = 0, p\ = 0 and p\ = ~ |i| 3 (pj,pj)} 
and is a submanifold of T*Jq. Hence, range 2zf is a submanifold of T*Jq . q.e.d. 


Recall that the Liouville form of T*J ] {) is 


6 = p t dt + p x dx + p^ dx 


(95) 


with exterior derivative 


o) ~ dh ~ dpt A dt + dp x A dx + dp x A dx (96) 

the canonical symplectic form of T* . 

For / € C“(rV‘), the Hamiltonian vector field of / is the unique vector field 
Xf on T*/q such that 

Xf—1 at = -df. 
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where _l denote the left interior product (contraction). The Poisson bracket of 
two functions /i,/ 2 e C co (T*J ] Q ) is given by 

{fuf2}=X A (fi). (97) 


It is bilinear, antisymmetric, and it satisfies the Jacobi identity 

For the sake of future convenience, define the reparametrization-invariant func¬ 
tion h by 


i W / v , <P- 

h = —r Px) + —n 

4 \x\ 


x) 


(98) 


Note that h is smooth, because x f 0. and we can use the constraint h = 0 instead 
of |x| 3 (px, px) + 4 (p x , x) = 0 in describing the range of Jz?. In other words, 


range = [it, x, x, p t , p x ,p x ) eT*J l 0 \p, = 0, (p x , x) = 0 and h - 0). (99) 


The Hamiltonian vector fields of the constraint functions p t , (p x , x) and li are 




Pt 


X (Px4) 


x h 


d_ 

dt 

d d 

x—-Px^~, 
OX OPx 


1 .. d 1 . <9 

-(x,x)p x -- + tt:X— 

2 ox \x\ ox 


1 9 ( 1 / 

~ T\P*^— + (P*’P*> + 

\x\ op x \ 2 


(Px, x) 

|ip 


d 

X dp.x 


Note that all the Poisson brackets of the constraint functions vanish identically 


{(Px,x),p t ) = i h,p t ) = {h,(p x ,x}} = 0. (100) 

This implies that range Jz? is a coisotropic submanifold of (T*jLa>). 


4.2. Action of Diff + R on T*J^. Recall that for X = t3, e diff + R, the action of 
the one-parameter subgroup exp sX on is generated by the vector field X 1 = 
t - tx-gr. The lifted action of exp sX on T*J ] {) is generated by the Hamiltonian 
vector field X j t , where 

C / T (0 *> •*> Pit Px, Px) = < Ptdt + Pxdx + p x dx, X T (t, X, x)> = T(t)p, - fit) {Px, x). 
The map 

d 

/m : diff+R -» T*j\ = T (t)p, - f(t)(p x ,x) 

may be interpreted as the momentum map for the action of the group Diff + R on 
T*Jq. Writing it this way avoids unnecessary discussion about the topology of the 
dual of the Lie algebra diff + R. The constraint equations p, — 0 and (p x , x) = 0 
imply that cy fditr vanishes on range 2. In other words, 

range ££ c /^(O). 

Proposition 4.3. ^ di( ](0) is a coisotropic submanifold ofT*J q. The null distribu¬ 
tion of the pullback of a> to #j7^.(0) is spanned by the Hamiltonian vector fields 
Xp t and X(p x xy 

Proof This follows from the proof of theorem l4~2l and equation (1100 1) . 


q.e.d. 






40 


L. BATES, R. CHHABRA AND J. SNIATYCKI 


Integral curves of the Hamiltonian vector field X Pt = arc lines parallel to the 
/-axis. Integral curves of X( Px j) satisfy equations 

d 

—x(s) = x(s), 

ds 

d 

TPx(s) = ~Px(s). 

ds 

Hence, for each p = (t,x,x,p t ,p x ,Px) £ T*/, 1 ,, the integral manifold of the distri¬ 
bution on T*Jq spanned by X Pt and X< Px ^) that passes through p is 

O p = {{u,x,e s x,pt,p x ,e~ s px) \ (u, s ) e R 2 }. (101) 

Theorem 4.4. For each p e c T*J p , the orbit of the Lie algebra diff + R 

through p and of the reparametrization group Diff + R coincides with the integral 
manifold O p given by equation (1707 D . where p, = 0. 

Proof Orbits of the action of the Lie algebra diff + R on T*j\ are orbits (acces¬ 
sible sets) of the family { Xje T r| 6 diff + R} of Hamiltonian vector fields on 
7*7,1. Since X = X Tp< - Xp Pi ^) and J? T vanishes on ^FgCO), it follows that the 
restriction of Xjr T to ,^^,[(0) is 

= rX P,\/ A fO) ~ fX iP*Mf A fOY 

Therefore, / T \/-. l ( 0 ) ' s contained in the distribution spanned by ^ Pl \ j 1 ( 0 ) an d 
X{ Pi x)\ r /-'( 0 y Hence, for each p e ^^,[.(0), the orbit of diff + R through p coincides 
with the integral manifold O p given by equation (11011) . 

The reparametrization group Diff + R acts on 7 ( * by 


Diff + Rx/g^ : (<p, ( t , x)) i-» 

where tp is a smooth function on R such that fit) > 0 for t e R The lift of this 
action to T* .7,' is 




Diff+ R x T*J' -» T*j\ : (<p,p) i-» 


4 p, ( Px,x)ip(t) 


It preserves #, if J.(0). Hence, the orbit of Diff + R is 

{{p{t),x,f~ x (t)x,f(t)~ x p t ,p x ,f(f)px) e T*Jq | <p e C°°(R), <p(t) > 0}. 

The action of Diff + R preserves #^(0), given by p t = 0 and (p x , x) = 0. Hence, 
orbits of Diff + R contained in ^,,](()) arc 

{((pit), x, (p~ x (t)x, 0, p x , (fi(t)px) € T*j' 0 | (px,x) = 0, p e C°°(R), <p(t) > 0). 


For each t, pit) = u and (pit) = -s, are independent. Therefore, orbits of Diff + R 
contained in #^(0) coincide with the corresponding integral manifolds given by 
equation (11011) . q.e.d. 
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4.3. Reduction of Diff + R symmetries. In this section, we discuss the space 

R = J, ^(0)/Diff + R 

of Diff + R-orbits in c / r ,T ( |(0). According to Theorem I4.4[ the reduced phase space 
R is the space of integral manifolds in #^(0) of the distribution spanned by X Pi 
and X( Px ^). We have shown that the orbit of the vector fields [Xp^X^.^f through 
p = (;t, x, x, p t , p x , px) € ^(0) is 

O p = i (u, x, e s x, p t , p x , e~ s px) | (; u, s) € R 3 }. (102) 

We are going to show R is a quotient manifold of ^7^(0), which will imply that 
R has a unique symplectic form ojr such that 

p*co R = Cut, (103) 

where l : is the inclusion map. 

In order to parametrize the reduced phase space R, define spherical coordinates 
(r, a, ft) by 

x 1 - rsinftcosa, (104) 

x~ = r sin ft sin a, 

x 3 = r cos ft, 

together with the dual momentum variables {p^, pp, /;„) defined by 

Pxdx = pt-dr + ppdft + pada. (105) 

Proposition 4.5. p,- = (p x , x) /r. 


Proof. This is a simple verification. 


q.e.d. 


Denote by p : ^,^(0) —» R the reduction map associating to each point in 

r /:; f (0) the orbit of \X Pi , X( Pr j )! through that point. Let 

5 = \(x, x) € TR 3 | |jfc| - 1} 


be the unit sphere bundle over R 3 parametrized by coordinates (x, a, ft). The Liou- 
ville form of T*S is 


and 


% = Pxdx + ppdft + Pada, 


(106) 


<4>s - d6s 

is the canonical symplectic form of T*S. 


Proposition 4.6. There is a unique symplectomorphism x : (R, ojr) —> (T*S,a>s) 
such that 

xo P : J. dTff(0) T*S : (t, x, x, 0, p x ,px) ^ (x,ft,a,p x ,pp,Pa), 

where the (ft, a, pp, pf) are related to (x, x, p x , Px) by equations U04\) and \105\t . 
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Proof. Consider first the space R\ = p t 1 (0 )/X Pt of integral curves of X Pt in p t 1 (0). 

It is a quotient manifold of p~ l ( 0) with projection map 

pi : pj\0) Ri : (t,x,x,0,p x ,p x ) (x,x, p x ,p x ). 

Moreover, it is a symplectic manifold with the symplectic form 

(i)\ = dp x A dx + dp x A dx. 

The constraint function (p x , x) is left invariant by the action X Pt , and pushes for¬ 
ward to a function on R\, denoted by (p x , x) ,. That is, (p x , x) = pj (p x , x) x . More¬ 
over, the Hamiltonian vector field X( PitX ) restricted to p~' (0) pushes forward to 
the Hamiltonian vector field on R i corresponding to the function (p x , i), on R\. 
Denote this vector field by X( Pi>x y l . 

By definition, r = \x\ + 0 on T*Jq. Since p x r = (p x ,x), it follows that on 
pi l/^m) the Hamiltonian vector field of (p x , x)\ is proportional to the Hamil¬ 
tonian vector field X [>r = Jr. Therefore, the space A '2 = {p x ,x)\ x (0)/X< /)i;i > 1 of 
orbits of X( Pi ' X ) l in (p x , x)~[ x (0) can be parametrized by (x,fi,a, p x , pp, pf). It is a 
symplectic manifold with the symplectic form 

u >2 = dp x A dx + dpp A dfe + dpa A da. 

The coordinates (x,0, a, p x , pp, p„) define a symplectomorphism between (Rj, a> 2 ) 
and (T*S, a>s), where u) s is the pullback to T S of the canonical symplectic form 
on r*(TR 3 ). However, R = )/{X Pt ,X( PuX )} with the symplectic form ojr 

is naturally symplectomorphic to (R 2 ,cl> 2 ). Hence, (R. ojr) is symplectomorphic to 
(T*S,a>s). q.e.d. 

It follows from Proposition 14. 61 that we may identify ( R , ojr) with ( T*S,a>s )■ 

The action of the Euclidean group SE(3) on R induces a Hamiltonian action 
of E on TVq generated by the Hamiltonian vector fields X Px , X P/j and X p . t . This 
action preserves the constraint functions p t , (p x ,x) and h. In particular, it induces 
an action of SE(3) on the zero level set of the momentum map for the 

action of diff + R. On the other hand, the action of E on R 3 induces a Hamiltonian 
action of SE(3) on T*S , presented as 

T*S - I (x,X,p x ,p x ) € r*(TR 3 ) I |x| = 1, Pr = 0), 

which is generated by the Hamiltonian vector fields of p x , pp, and p (t considered 
as functions on T*S. Moreover, these actions of SE(3) are intertwined by the 
reduction map p : —> R followed by the identification R = T*S. 

4.4. Hamiltonian dynamics. The range of the Legendre transformation is char¬ 
acterized as 

range «£? = 0 ) n /7 _1 (0), 

where 

IM 2 , . , (Px, X) 
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The Hamiltonian vector field of h is 


v 1 /• -x 5 1.5 1 

X h = -(x,x)p jc — + - - -A- - - - — 

2 dx (x,x) 112 dx (a, i} 1 ' 2 


d ( 1. 

Px ~dp \2 ( ' Pi,P ^ + 


(p x ,x) \ . d 
a.v>- 5 2 / v 5/t' 


In order to find the integral curves of X/, on the range of «£?, observe that they 
satisfy the equations 


d 





(x, x) 


1/2 


X, 


- (x, X) Px, 


( X, x) 


1/2 


Px + 


1 / \ , (Px,x) 

-2 <p, ’ p ‘ ) + dJy 371 


x, 


0, 

0 . 


Multiplying equation j^x = \ (x, x) p x by x, and using the constraint equation 
( px, x) = 0, yields ^ (x, x) = 0, hence (Ms), Ms)) = |io| 2 . Since h is reparametriza- 
tion invariant, without loss of generality, we may assume that |io| = 1, that is, the 
arclength parametrization. Moreover, on the range of Jzf, li = 0, which implies 
that \ (px, p^ = -2 (p x , x). This leads to 


li 

x. 

(107) 

d . 

1 

(108) 

Ts x 


d 

d~s p * = 

-Px + 3 (p x , x) X, 

(109) 

ii 

£ 

o, 

(HO) 

d 

ds 

0. 

(111) 


Equation (II 10b implies that p x is constant. The angular momentum is given by 
l - x X p x + x X px (see equation ([6T1 ) Hence, 




x p x + xx 




x p^ + x x 



x x Px + -Px XPx + XX (-p x + 3 (Px, x) X) 

0, 
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and so is conserved. Therefore, (p x , l > = (p x , x x p x ) is also conserved. Multiply¬ 
ing equations (1 1081) and (1109b by p x yields 

d 

— (p x , x) = (p x , p x ), 
ds 


T (Px, Px) = - (Px, Px) + 3 (p x , x ) 1 , 

ds 


or 


(Px, x) = - (p x , p x ) + 3 (p x , x) 2 . 


Multiplying by j^(p x ,x) and integrating gives 


1 / d < p x , x) 

2 \ ds 


= ~ (p x , Px) (Px, X) + {p x , i> 3 + constant, 


which can be integrated since it is separable. If p x d 0, then this equation gives the 
component 

(Px, X) 


X" = 


I P: 


2 'Px 


of x parallel to p x . Integrating x^(s) yields the component of the motion in the 
direction of p x . Returning to equations (1108b and (1109b gives 


Ts X = 
d 

-rPx = ~Px + 3 {p x , x) x, 
ds 

where (p x , x) is assumed known from the discussion above. Hence, 

d 2 . 1 3 

tf x = -2 p ‘ + 2 Mx - 

Writing x and p x in terms of their components x' and p x ,, 

d 2 1 3 .. . i 

= -~p x i + -(p x ,x)x. 


Division by x\ implies 


^ I in 

ds 


ds 


1 d 2 x‘ p..i 3 , 

= + n<P x ’X>> 


d£ ds 2 2— 2 


ds 


ds 


which implies 


In 


or 


ds 

In 


■If 


(p x , x) (s) ds 


xPxi/ 2 — x‘ 

ds 


(p x , x) (s) ds, 


( 112 ) 

(113) 

(114) 

(115) 


(x‘) p x ,2 —x l - cexp( | f(s)ds |, 


f 


so that 




















ELASTICA AS A DYNAMICAL SYSTEM 


45 


where c is a constant dependent on the initial data. Integrating once more yields 


1 


c + 1 


CO 


i\t+p xi /2 


if p x i / 2 -1, and 

if p x i 12 = -1. 


In lx I = 


= cjexp(j 

= J“ P (/ 


/(ds 


/(s)ds ds 


Remark 4.7. The Hamiltonian vector field of p t + h is 


X, 


d 1 d 


1 d 


Pi+h 


- Xp, - Pl - + jTp - jTp + 3 


. d 

|d| 5 X dp x 


In the arclength parametrization, it is 


ds 

1, 

ll 

X 

~ ^ l"Q 

X, 

d . 

1 

Ts x 

2 P *’ 

d 

d~s Pt = 

-Px + 3 (p x , X) X. 

d 

0. 

-rPx = 

ds 


Hence, t = to+s, and the solutions of this system can be obtained from the solutions 
for integral curves of X/, by replacing s by t - to- 


Theorem 4.8. If (to, xo,xo,Q, p X0 ,Px 0 ) - -£f(to, *o, *o, xo, xo)> the solution of the 
Euler-Lagrange equation with initial data (to, xo, xo, xo, *o) is equivalent to the in¬ 
tegral curve ofX Pt+ k through (to, xq, xq, 0, p X0 ,Px 0 ) given above. 


Proof The Euler-Lagrange equations for a second order Lagrangian are 

dL d dL + d 2 dL 
dx dt dx dt 2 dx 


Since p k = §§ and p* = § k- f t p k , they are §| - j,Px = 0. For elastica. 


L(x, x, x ) 


|.x| 2 (x, x) 2 

W 1 " W 5 


H = 0, and the Euler-Lagrange equations reduce to 

d 

7P = °- 

Since p x = - j t p x , the evolution of p k is given by 

d dL 
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where 


Px = 


dL „ x „ (x, x) x 

_— 2_2—_-_ 

dx |jfc| 3 |x| 5 


and 

dL _ d !\x\ 2 {x,x) 2 \_ \x\ 2 . 0 (x,x) .. (x, x) 2 , 

<9x <9x \ |x| 3 |x| 5 / |x| 5 |x| 5 |x| 7 

In the arclength parametrization \x\ — 1, (x, x) = 0, px — 2x and 


dL 

dx 


= -3 |x| 2 x = -- Ip.vl 2 x = 3 (p x , x) x 


because 

h = - (Px, Px) + ( Px , x> = 0 

on the range of _§?. Thus, the Euler-Lagrange equations of elastica in the arclength 
parametrization are equivalent to 


Jt Px 


dt Pi 


—x 

dt 

d 

—x 

dt 


0 , 


Px + 3 (p x , x) X, 


1 

2 Pi ’ 


X. 


This system of equations, together with the substitution t = to + s, leads to the 
equation for integral curves of the Hamiltonian vector field of X Pt+ h given in the 
remark above. q.e.d. 


5. Quantization 

There are two ways of quantization of a system with constraints: the Bleuler- 
Gupta quantization of the extended phase space followed by the quantum reduction 
liTI . |8| , and Dirac’s classical reduction followed by quantization of the reduced 
phase space studied by Dirac in |4). Here, we follow the Bleuler-Gupta approach. 

The extended phase space T*J ] () is the cotangent bundle of an open subset of R 14 . 
Therefore, it is convenient to use the geometric quantization in terms of the vertical 
polarization tangent to fibres of the cotangent bundle projection. This approach 
leads to Schrodinger quantization. Since the functions we want to quantize are at 
most quadratic in momenta, we can use results on Schrodinger quantization as they 
are presented in texts on quantum mechanics. For technical details see fl9l . 

Consider the space C,“(./ ( I ) )8C of complex-valued compactly supported smooth 
functions T on endowed with the scalar product 

(T't | TS) = j Thf/,x,xjT^ft,x, x)dt d 3 x d 3 x. (116) 

Jjl 

The completion of C^°(7 q)®C with respect to the norm given by this scalar product 
gives rise to the Hilbert space §o °f quantum states of the system. 
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If / e C°°(/q), then the operator Q n *f corresponding to the pullback n*f of 
/ by the cotangent bundle projection n : T*J^ —> acts on C°°(J ( *) ® C via 
multiplication by / 

QrfV = f'V. (117) 

Operators coiTesponding to canonical momenta act by derivations with respect to 
the coiTesponding variables in j\. In particular, the linear momentum p x gives rise 
to 

<ii8) 

and the angular momentum / = x x p x + x x p x leads to 


Qj = -itix x- ihx x —, 

ox ox 

where fi is Planck’s constant. Furthermore, ,/ T = T(l)p, 


(119) 

- i(t) (px,x) quantizes to 




+ ifif 




- itif. 


( 120 ) 


Moreover, for h = ^-j- (px,Px) + the metric on R 3 giving the quadratic term 
has scalar curvature 2, and the coiTesponding operator is 



iti I d \ 


( 121 ) 


where 

a -(s4) 

is the Laplace operator in the variables i'Q These differential operators on C°°(/q)® 
C extend to self-adjoint operators on §o- 


5.1. Quantization representation of SE(3). The skew-adjoint operators -pQ p . 
and -rQi generate a unitary representation of the Euclidean group SE(3) on $o 
such that for g e SE(3) and T 1 e C“(/q) ® C 

Ug'V (t,x,x) = 0*_ 1 l P (t.jr.jc), 

where 

O : SE(3) x 7’ ^ 

is the action of SE(3) on 7 ( '. The invariant vectors of this representation arc eigen¬ 
vectors of -jQ p . and ~pQi coiTesponding to the eigenvalue 0. In other words, in¬ 
variant vectors T 1 of U are characterized by the equations 

Q P v = 

Qi'v = o. 

Since SE(3) is not compact, invariant vectors of U are distributions on J ] () . 


4 Correction terms in equations 11201 and 11211 are consequence of using half-forms in order to 
ensure that the space §o of quantum states can be described as the space of square-integrable complex 
functions on yj (see Cl-) 
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5.2. Quantization representation of Diff + R. The reparametrization group Diff + R 
acts on C°°(7q) 8 C by the pullback of its action on J ] () 

Diff + R x (C“ (/,',) ® C) -> C 0 ^./, 1 ,) ® C : (<p, T) i-» ( y r 1 ) 1 *'F, 

where ip~ ] is the inverse of ip, and 

(<^r 1 ) 1 *'f'(t,x,x) = T fp~ l {t,x,x)) = x ¥\(p~ 1 (t),x, -—) . 

\ (<p )V)I 

For an infinitesimal diffeomorphism ip e (t) - t + er{t) + ... generated by T{t)d t , 


(^> e 1 ) 1 * v P(f, x, x ) = v F(t, x, x) — £ 




x ¥ + ... 


and 


—(ip e 1 ) u '¥(t,x, x)| e =o 
de 




v P(t, x, x) 


( 122 ) 



vj (t, x, x) - f'Vit, x, x). 


Therefore, 

(jQ a fr' ¥ )^X,X) 


-rUPe 1 ) 1 * X k(F -X, -f)le=0 + t'Ffc X, x) 
de 

^-{(^ 1 ) 1 * x P(t,x,x) + <j> £ (ty T(f, x, x ))}| e=0 
de 

-7- 1 ) 1 * [^e(0W, *)] }le=0 • 


This establishes 


Proposition 5.1. The operator 


—Q a = T— - iflT ( X, — ) + f 


dt 


dx 


generates the action on ® C of the one-parameter group <p e — exp er given 

by 

where 


H : (p e , , F) ^ 


x, x) = (^ e 1 ) 1 *[^ e 'P(t,x,x)]. 


Theorem 5.2. The map 

S : Diff+ R x (C 00 ^) ® C) -» C 00 ^) ® C : (<p, T') i-» E^T, 

where 

S x, x) = Jc, x)], (123) 

/.v n linear representation of Diff + R on C‘ X ( J ] {] ) 8 C preserving the scalar product 
given by equation (II 161) . 
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Proof. Clearly, E^ acts linearly on C°°( 7 q) 0 C. For ip\,p>2 e Diff+R and 'F e 
C°°(Jl) 0 C, 

- (^ i ) i *!^ 2 (0(^r 1 ) 1 *[^i'F]} 

- (^ 1 )>2(^2 1 )*(^r 1 )>i[(^2 1 ) 1 *(^r 1 ) 1 *'J'] 

- (^2 1 )*^2(^r 1 )>i][((^2 o ^i)- 1 )^] 


UfPi 0 <pi) l T(<pi ° ^i)(0 


(ip2°ip\){{ip2 0(piY l {t)) 

(p 2 ° ^ixc^r 10 v , 2 1 )W) 

(p2 ° <pix(.<pi\<p?m 

f>2 (<fi 1 (pi 1 (pi 1 (0)))<^ 1 c^r 11 (^r 1 (0)) 

v?2 (^2 1 co)^ 1 (^>r 1 (^2 1 co)) 


as required. 

For 'Fj, ^2 € C°°(/i) 0 C and (p 6 Diff+ R, 

(S^^i | E^ X F 2 ) = I E v '¥i(t,x,x)E ip '¥ 2 (t,x,x)dtd 3 xd 3 x 


I (ip 1 ) 1 *{^)(t) x Pi(t, x, x)} (jp l ) 1 *{ip{t)'¥^(t,x,x)}dtd 3 xd 3 x 

JP 0 

J" [(pitp-'m 2 x ¥i((p~ 1 (t), x, x/((p- l )(t )) {'Fi( y T 1 (/), x, xl{if~ l )(t)} dtd 3 xd 3 x. 


Note that the inverse function theorem guarantees 

<K<P~' 1 ( 0 ) = 7—77—- 

(p~ l m 

Introducing new variables 

t = p>~ l {t), x — x and x - x/(p~ ] )(t), 

yields 

dt = (p^ l )(t)dt, 


dx = d.x, 
dx = 


1 (<f l )(t)x u 

dx -—-r dt , 


1 


((9 'tof 


■ dt d 3 xd 3 x. 


dt d 3 xd 3 x’ 
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so that 

Therefore, 

(H^l I 5^2) 


dtd 3 xd 3 x = [(</> ')(?)] dtd 3 xd 3 x. 

j [ip(tp~ 1 (t))] 2 , ¥ 1 (t, X, x') x I J 2(t, x, x') dt d 3 xd 3 x 

j -:- x, xfy¥2(t, x, x') [(^> -1 )(0l dtd 3 xd 3 x 

[(^ _1 )(0] 2 L J 

I T| (?, x, x') l Vi(t, x, x) dtd 3 xd 3 x. 

Jj‘ 

q.e.d. 


Note that, for T, ¥' e C“C/q) 0C, the integral defining the scalar product (II161) . 
and 

(T' | T) = f W'(t, x, xfY(t, x, x) dtd 3 xd 3 x (124) 

J j b 

can be interpreted as the evaluation on T of the generalized function (distribution) 
T" e (C,“(/(',) 0 C)'. The representation 3 of Diff + R on C ( “ (■/,') 0 C extends to a 
representation of Diff + R on (C 0o (7 ( 1 ) ) 0 C)', which we also denote by 3, such 

(H^T" | TO = Of" | 3^-.¥) (125) 

for ip € Diff + R, T" € (C™(J 3 ) 0 C)'. Note that, if T' € C°°(J l 0 ) 0 C, then the 
definition of the action 5^ on T", given here, coincides with the definition given in 
equation (11231) . 

It remains to examine the space of Diff + R-invariant functions. Since the group 
Diff + R is not compact, the only compactly supported Diff + R-invariant function in 
C“(7 q) 0C is identically zero. Hence, Diff + R-invariant functions are in C“ J (./,',) 0 
C. 


Lemma 5.3. For q = (to, xg, i'o) € jL the orbit exp(diff + R)(<jr) o/diff + R through q 
coincides with the orbit Diff + R( q) o/Diff + R through q 


exp(diff + R)(<jr) = Diff + R(q). 


Proof. For X = rd t e diff + R, the integral curves of X 1 = tJj - txjr satisfy the 
differential equations 

dt dx dx 

—— — r, — = 0, and — = -ix. 

ds ds ds 


Hence, 


and 




dt' 

T(t') 


= S. 
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Choosing r(t) = e 


f'—= (V*=,- 

Jtn T ( f ) Jtn 


e‘°. 


Hence, e‘ - e to = s, which implies that 

t = log |s + e~ ,0 \. 

But the range of the logarithm is (- 00 , 00 ). Therefore, the range of values of t on 
the orbit of dilf + R through q is (- 00 , 00 ). 

Moreover, for i = 1,2,3, 

dxi 

— = ~TXi 
as 

implies 

dxi 

— = -f(t(s))ds 

X{ 

so that 

x(s) = X( } exp i-r t(t(s)) ds 


Since r(t) is an arbitrary function of t, it follows that the orbit of diff + R through 
q = (t 0 , X{), xq) is 

exp(diff + R)(q) = {(u,xo,e v x) | (; u,v ) e R 2 }. 

The action of Dilf + R on is 

Diff+ Rx/j^ Jq : (<p, (t, x, x)) i-> ^ip(t), x, , 

where ip(t) > 0. Since, ip(t) and Cp(t) are independent, it follows that the orbit of 
Diff + R through q is 

Diff + R(g) = {(U, Xq, wx) I (u, w) £ R 2 , w > 0). 

Hence, exp(diff + R)(g) = Dilf + R(g). q.e.d. 

Theorem 5.4. A function T" e C°°(J^) 8 C is Diff + R -invariant if and only if 


Qjy = 0 


for all t e diff + R. 


Proof If T" e C°°(/q) 8 C is Diff + R-invariant, then it is invariant under the action 
of every one-parameter subgroup of Diff + R. By proposition 15.11 actions of one- 
parameter subgroups of Diff + R on C°°(Jq) 8 C arc generated by yQ,/ t f° r T € 
diff + R. Hence, Qjff" = 0 for all r e diff + R. 

Conversely, suppose that T' is a function in C°°(Jq) 8 C such that QjfV = 0 
for all t e diff + R. Hence, T' = T' for every one-parameter subgroup ip £ of 
Ditf + R. Recall that, for ip e diff + R, 

t,x,x) = (ip~ l ) u [<f(tY¥'(t,x,x)\ = [(ip~ l )*<f(t)Mp~ X ) U ^\t, x, .*)] 

= ip(ip~ ] (t))Y(ip~ ] (t, x, xj) = Y(ip~ x (t, x, x)), 

dt 
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because 

implies 

Therefore, 


p o p 1 = identity 


-l ,^ d P l (t) 

<p(<p (0)—— = 1. 
at 

S -iT J '(t,x,x) = -i-'P '(tp(t,x,x)). 


V- <p(t) 

It follows from the lemma above that there exists a finite sequence of one- 
parameter subgroups ip ei ,..., p €k such that 

pit, X, X) = p ek {. . . (p £l (t, X, X)) ...). 

Since 


and 


Hence, 


p(t,x,x) = \ p(t),x, —j , 

<f(t) = <Pe k i- ■ ■ ...) = Pe k ° ■ ■ ■ o p ei (t) 


= <p £k (. . . {p £l (t)) p ei (t) = ~^Pe k ° ° ^£,(0- 


H -1 Y(t, x, x) = -i- v F , (^(t, x, x)) = -y-'P' [pit), x, -p- 

V (Pit) Pit) 


1 


-r 


<£(0 


Pe k °---° <f ei it), X, - ---— 

Jt p €k o • • • o p ei (t ), 


° • • • ° MO 

= x, x) = E^ 0 ,„ 0 ^Y(t, X, x) = 3 ¥) _ ei0 ... 0 ^'P , (t, X, x) 

= 3 v _ e| • • • H^'P'O, X, x) = Y it, X, x), 

because T" is invariant under the action of onc-paramctcr subgroups of Diff + R. 

q.e.d. 


5.3. Quantum implementation of constraints. In the Bleuler-Gupta approach, 
the extended phase space is quantized first. This associates to each constraint func¬ 
tion the corresponding quantum operator. The next step is to implement the con¬ 
straint conditions on the quantum level. This is done by placing a restriction on the 
states of the system. 

Definition 5.5. Admissible quantum states are eigenstates of the quantum operators 
associated to the constraint functions corresponding to the joint eigenvalue zero. 

For elastica, the classical constraints arc 

Pt = 0 , 

(Px, x> = 0, 

|x| 3 (Px, Px) + 4 (Px, x> = o. 
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However, linear combinations of these functions with smooth coefficients leads to 
further functions that vanish on the range of Legendre transformation. We might 
not be able to quantize these functions in our chosen quantization scheme. This is 
why it is important to have criteria that help select a convenient basis of the ideal 
of functions that vanish on the range of Legendre transformation. 

The first two constraint conditions p t = 0 and (p k , x) - 0 arc equivalent to van¬ 
ishing of the momenta :: / T for the action of one-parameter subgroups of Diff + R. 
Therefore, the quantum implementation of these constraints is the requirement that 
the admissible wave functions T should satisfy the conditions 


Q = 0 for t e diff + R. 

By the results of the preceding section, this is equivalent to requiring that admissi¬ 
ble states should be invariant under the action of the quantization representation of 
0 C. 

The third constraint function has no immediately clear geometric interpretation. 
We have used the freedom of the choice of generators of the ideal of constraint 
functions and replaced it by the reparametrization invariant function 

, |d| 2 . (p x , x) 

n = — (Px, Pi-) h -. 

4 yix,Fx/ |i| 


The quantum operator corresponding to li is 


tr 


Q h = ~ 


i 


ml . a\ 

\M \ dx)' 


Thus, admissible states T of quantum elastica arc also required to satisfy the equa¬ 
tion Q^V = 0. 


Summary 5.6. The admissible states of quantum elastica arc given by functions 
T e C°°(/q) 0 C that arc invariant under the quantization representation of the 
reparametrization group Diff + R and satisfy the equation 




A _I 

3 


T 


®/,,AW = 


1*1 X ’ dx I 


'V = 0. 


As we have mentioned before, admissible functions of quantum elastica arc not 
square integrable on J ] {] . Therefore, in this formulation of the theory, we need to 
introduce a new scalar product on the space of admissible states using physical or 
geometric criteria. For example, we could use Diff + R invariance of admissible 
quantum states and relate them to smooth functions on the space S of Diff + R - 
orbits in , which satisfy the differential equation obtained from the quantization 
of h considered as a function on T S. 
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